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In this 

study , 

the 

Ground 

Hydrologic Model 

( GHM) 

developed by 

Lin , 

Alfano 

and Bock (1978) for use 

in an 

atmospheric 

general 

circulation 

model (GCM) has 

been 


refined. A series of sensitivity studies of the new version 
of the GHM were conducted for the purpose of understanding 
the role played by various physical parameters in the GHM. 

This version of the GHM has made the following 
refinements : 

1 ) The GHM is coupled directly with the planetary 

boundary layer (PBL), using Deardorff's (1972) 

parameterization . 

2) A bulk vegetation layer is added with a more realistic 
large-scale parameterization. 

3) The infiltration rate is modified using Green-Arapt's 


( 1 91 1 ) formula . 
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The Gh'M has been tested using input data derived from a 
GCM simulation run for eight North America regions for 45 
days. The results are compared with those of the resident 
GHM in the GCM. The daily average of grid surface 
temperatures from both models agree reasonably well in phase 
and magnitude. However, large difference exists in one or 
two regions on some days. The daily average 
evaportranspiration is in general 1 0-30^ less than the 
corresponding value given by the resident GHM. 

Sensitivity studies have been Conducted for: 

1 ) I/iitial .conditions and lower boundary conditions. The 
effect of different initial soil moisture conditions in the 
surface layer persists approximately one week, while in the 
lower layer at least as long as the operational period. 
Different lower boundary conditions only produce minor 
effects. 


2) Vegetation density. For extreme cases of 
desertification and afforestation, the effects of the 
vegetation density on surface temperature, evaporation, 
sensible heat and soil moisture are most significant. 
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3) Canopy resistance. Mainly the daily averaged grid 
evapotranspiration and the soil moisture content are 
effected . 

4) Surface albedo. Albedo exerts a large effect on the 
energy balance and thus temperature variation. 

5) Depth of root zone and root density distribution. The 
effects from these two parameters on the grid evaporation, 
sensible heat, moisture content and temperature are small. 
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heat conductivity of sand and silt (cal/cm/s/°K) 

heat conductivity of clay (cal/cm/s/°K) 

heat conductivity of organic matter (cal/cm/s/°K) 

air density (g/cm ) 

water density (g/cm^) 

bare soil fraction in grid 

variable in Equation ( 3 . 5 . 19b) (°K) 

vegetation fraction in grid 

soil matrix potential (cm) 

soil matrix potential in Table 2-4-1 (cm) 

86400 (sec) 


Note: All variables preceded by 'grid' indicate these 
variables are grid averaged. All variable preceded by 
•daily' indicate these variables are daily averaged. 
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Chapter I 
INTRODUCTION 


1.1 SIGNIFICANCE OF LAND SURFACE PROCESSES 

In the earth environment, there are two basic transfer 
processes, the transfer of water and heat between land sur- 
r . faces and the atmosphere. The hydrologic cycle deals with 

the transport of water in the gas, liquid and solid states. 
Precipitation as one of the main components in the cycle is 
the predominant source of water supply to the earth. It 
contributes to >. . soil moisture and groundwater through in- 
filtration and deep percolation, or remains as surface wa- 
/ 

ter through surface or subsurface runoff. However, approxi- 
mately 60 to 70# of precipitated water on land surfaces 
returns to the atmosphere through evapotranspiration. The 
part of the hydrologic cycle near the atmosphere-land in- 
terface can be isolated and described by the principles of 
conservation of water mass. For a one-dimensional soil 
column of finite depth below the ground surface, the equa- 
tion for the water balance is 

dM/dt=P r -E-R Q -R b (1.1.1) 

where M is the water stored in the column, P r the precipi- 
tation including snow and ice melt on the land surface, E 


1 



2 


the evapotranspiration, R Q the surface and subsurface ru- 
noff, and R^ the deep percolation. All variables in Eq. 
(1.1.1) except P r are closely related to the soil moisture 
content in the column. 

The other important process on land surfaces is the heat 
exchange between land surfaces and the atmosphere. Short- 
wave solar irradiation and long wave irradiation from the 
atmosphere are the main supplies of the thermal energy for 
the hydrologic system. The latent heat plays an important 
role in linking the thermal energy exchange with the water 
cycle . 

If the heat transfer of snow and ice melt and plant pho- 
tosynthesis are neglected, the net heat flux to the ground 

G can be determined by the energy balance on the ground 

/ 

surface , 

G=R n - L*E - H (1.1.2) 

where R n is the net radiation, L the latent heat of vapori- 
zation of water, E the evapotranspiration upward from the 
surface and H the sensible heat flux. The terras • in the 
right-hand-side of Eq. (1.1.2) are greatly affected by 
ground surface temperature and in turn the surface tempera- 
ture is determined by the heat transport process. 

The total evaporation in a year averaged over the land 
is 0.45 ro (Brutsaert, 1 981 ) per unit surface area. It is 
much more than the fresh water stored in a soil column 



(O.lftl-m) or in an atmospheric column (0.0274 m) of the 
same cross-sectional area. In other words, the hydrologi- 
cal process at the land surface is quite capable of renew- 
ing the moisture in soil and in the air. 

In recent years, it has been shown in the course of at- 
mospheric general circulation modeling that land surface 
processes exert a great influence on atmospheric movements 
and climatic changes (Eagleson, 1 981 ; Mintz, 1981; Shukla 
and Mintz, 1982). Numerical experiments have been carried 
out under different land surface conditions and parameteri- 
zations. The results show large differences in the calcu- 
lated precipitation, temperature and motion of the atmos- 
phere. These studies have indicated that a physically 
realistic land surface hydrology model is essential to pro- 
vide a, GCM with a better estimate of exchange of pertinent 
properties on the atmosphere-land surface. 

Mathematically, a GCM needs to be supplied with lower 
boundary conditions at the land surface. An interactive 
ground hydrology model should be designed to simulate the 
soil moisture content and ground temperature as well as 
estimate the moisture and heat flux6s across the atmos- 
phere-land interface for the purpose of providing the 
source terms in the governing equations of GCM. 

GCMs have been developed through several generations 
(Carson, 1981) and have now reached a relatively advanced 
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stage. The most current GCMs are complex an<T make increas- 
ingly large demands on the most advanced computers. How- 
ever, their representation of land surface processes is 
still crude and simplistic. This situation urges research- 
ers to engage a more active and concerted effort on the 
land surface process study in order to provide a physically 
realistic GHM for use in GCMs. 

1 .2 PREVIOUS STUDY 

Reviewing the previous work, one finds that the study of 
land surface processes interactive with a GCM has gone 
through a number of developmental stages. In the earlier 
stage, GCMs use only fixed boundary conditions, for exam- 
ple, completely saturated soil with zero heat capacity. In 

fact, there was no ground hydrology involved. The concept 

/ 

of zero soil heat capacity was used by Manabe et al. (1965) 
and later employed in a GCM by Gates and Schlesinger 
(1977). The consequence of this assumption is an enlarge- 
ment of the amplitude of diurnal surface temperature varia- 
tion. A 'bucket* hydrologic model was developed by Manabe 
(1969) and used by Washington and Willimson (1977). This 
model, was designed to incorporate the most critical feature 
of soil moisture without resorting to a more rigorous ap- 
proach such as Philip and deVries (1957). The depth of wa- 
ter contained in the bucket represents the soil moisture 
per unit surface area in the plant root zone available for 
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evupotranspi ration . The bucket has a capacity equivalent to 
the maximum available soil moisture in the root zone or the 
field capacity of soil. By following Budyko (1961), actual 
evaporation is scaled from the potential evaporation 
through a proportionality factor which depends on the 
available soil moisture content. Runoff occurs only if wa- 
ter contained in the bucket exceeds the depth of the buck- 
et. In the model, the effects of different land surface 
cover, different soil type, and their spatial and temporal 
variation are not considered. Thus, the bucket model ap- 
pears to be oversimplified. Deardorff (1978) proposed a 
hydrologic model with inclusion of a canopy layer. In his 
model, the force-restore method was used to describe soil 
temperature, which is distinguished by a ground surface 

temperature with diurnal variation and a layer-average 

/ 

temperature with seasonal change. Vegetation canopy were 
modeled as a single layer with no heat capacity. Soil mois- 
ture was calculated in a way similar to the ground tempera- 
ture. The soil is divided into two layers: a upper surface 
layer of 10 cm thick and deep layer of 50 cm thick. The 
equations describing moistures in these two layers are sim- 
ilar to those describing the temperatures but with differ- 
ent physical parameters. 

The model developed by Lin, Alfano and Bock (1978) is 
one of the more advanced models for use in a COM, which in- 
cludes physically-based hydrological phenomena. Since this 
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study is developed from this model, it will be discussed in 
some detail. 

It is a two-layer ground hydrology model, which was de- 
signed to dynamically interact with the atmospheric general 
circulation model developed at the NASA Goddard Laboratory 
for Atmospheric Science (GLAS) as described in Halem et al. 
(1979)» Somerville et al.-(1974) and Tsang and Karn (1973)* 
Figure 1.2.1 (from Alfano, 1981) demonstrates schematically 
the main physical features of the model. 

In order to be compatible with the GLAS GCM, the land 
surface grid formation and the grid corresponding code in 
this model are exactly the same as the GCM (see Tsang and 
Karn, 1973)* Each grid is a rectangular cell of 4° longi- 
tude by 5° latitude. -The grid code for the latitude and 
longiWde is shown in Table 1.2.1. 

A global characterization of soil and vegetation was in- 
corporated into the model. Five types of soil texture from 
fine sand to clay loam based on the work of Ritjeraa (1970) 
and twelve categories of' vegetation from tropical rain for- 
est to desert (Strahler, 1971) were used. A soil texture 
and a' vegetation type were specified for ..each grid with 
further subgrid parameterization of vegetation density by 
fractional portions of bare soil and vegetative cover. 
Their relative fraction changes with grids and is also ad- 
justed for seasonal variation and for latitude. In this 
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model, the ground temperature was calculated by the modi- 
fied force-restore method (Lin, 1980). Even though the root 
zone was layered similar to Deardorff'a, the rate equations 
for soil moisture in the two layers were based on the Rich- 
ards law, in which all coefficients are derived from physi- 
cal principles. 

The prognostic variables for describing the soil mois- 
ture in this model are volumetric moisture contents in the 
surface layer and in the lower layer and for evaluating the 
heat budget are grid surface temperature and layer average 
temperature. The specification of the two layers for soil 
moisture movement is not necessarily the same as those for 
heat transport. The outputs from the model also include 
latent heat and sensible heat fluxes. 

Recently, a multi-layer canopy model was proposed by 
Mintz et al. (1983) for the NASA GCM, which refines to some 
extent the Lin et al. model to account for the effect of 
rainfall interception. 
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Table 1.2.1 Relating between (1,0) and (Longitude, Latitude) coordinate systems. 
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Figure 1.2.1 The ground hydrology model formulation 
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1.3 SCOPE OF THE STUDY 

In this study, several improvements on the Lin et al. 
model were implemented and a series of sensitivity studies 
of the GHM for a number of the pertinent physical parame- 
ters has been conducted. 

This new version of the model is designed to couple with 
atmospheric variables in the PBL in stead of the surface 
values which are derived by extrapolation from the results 
in the lower layers of the GCM (Alfano, 1981). This ap- 
proach will eliminate the crude approximation by the extra- 
polation. Moreover, in view of the large horizontal scale 
used by the GCM and GHM, it is inconceivable that the at- 
mospheric conditions can be uniform over a grid near the 
land surface, say, at the ’anemometer height’. The actual 
physical system has the PBL between the free atmospheric 
field described by the GCM and the land surface (Brutsart, 
1981). It is this layer that plays a dominant role in the 
transfer of momentum, moisture and sensible heat through 
turbulent transport and free convection. 

Although many schemes for evaluating the turbulent 
transport in PBL have been proposed, one can classify them 
into two categories.' Mellor and Yamada (1974) proposed a 
method using several layers within the lowest 2-3 above 
the land surface to resolve the vertical structure of the 
PBL explicitly. This method is clear and direct but it re- 
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quires several finite-difference levels in the PBL to 
achieve tolerable computational accuracy. 

Another approach is to parameterize all aspects of the 
PBL for a GCM. The drag coefficient and the bulk heat 
transfer coefficient are linked with a bulk Richardson Num- 
ber defined for the entire PBL. The accuracy of this ap- 
proach depends on the information available to describe the 
PBL. This knowledge has been accumulated recently. The pa- 
rameterization suggested by Deardorff (1972) has been test- 
ed numerically and appears to be adequate for use in exist- 
ing GCMs. The parameterization scheme will be adapted in 
this study. 

A bulk canopy layer over the vegetated ground surface is 

added to the model. Vegetation canopy plays an important 

/ 

role in the heat balance and water cycle and exerts signif- 
icant physiological and morphological effects on land sur- 
face processes (Deardorff, 1978; Mintz et al., 1983)* Fig- 
ure 1.3*1 (Mintz et al., 1983) describes the mean annual 
water and energy balance of a short and a tall vegetation, 
as derived from the measurements in two adjacent catchments 
having, nearly the same atmospheric conditions. The total 
evapotranspiration loss from the forest catchment is more 
than twice as large as that from the grass covered catch- 
ment. With the grass cover, 58£ of the net radiation energy 
on the surface is used for evapotranspiration and 42j£ for 



the sensible heating of the atmospheric boundary layer. But 
with the forest cover, the energy used for evapotranspi ra- 
tion exceeds the radiational heating by 15 %, and this re- 
sults in a removal of sensible heat from the boundary lay- 
er, i.e. a negative Bowen ratio. Although the atmospheric 
conditions are about the same for the two catchments, the 
net radiation R n is 10 % larger with the forest cover than 
with the grass cover. This is not only because the forest 
is darker and absorbs more of the incident solar beam, but 
also because the lower forest temperature will cause the 
forest to emit less infrared radiation. This example dem- 
onstrates that a canopy either modeled as a bulk layer or 
multi-layers between the atmosphere and ground surface are 
critical for describing land surface processes. 

The role played by infiltration in the partition of pre- 

/ 

cipitation is very important. Its rate can affect the mag- 
nitude of surface and subsurface runoff as well as the 
amount of water stored in the soil. The infiltration rate 
is determined by many factors such as vegetation type, soil 
texture, rain intensity. and so on. Its accurate evaluation 
requires sophisticate numerical methods. 

There are several methods being used currently to esti- 
mate infiltration rate. The most known are the empirical 
equation of Kostiakov (Childs, 1969) and Horton (1940) 
which have been popular because of their simplicity and ap- 
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plicabi.lity to most situations. However, both equations 
contain parameters that are difficult to predict because 
they have no physical basis. A more recent empirical equa- 
tion was given by Holtan (1961 ) . It expresses the infil- 
tration rate as a function of unoccupied pore space in 
soil. The model is convenient to use but its parameters in- 
cluding the characteristic depth are difficult to deter- 
rai ne . 

A simple equation proposed by Green-Ampt (1911) has re- 
cently attracted new attention because it has clear physi- 
cal meaning and can be expressed by a simple formula de- 
rived from the Richards law. In this study, this method 
will be adopted. 

Finally, the atmospheric stability effects will be in- 

/ 

eluded'. When the atmosphere is stable, turbulent mixing is 
suppressed and exchanges between the air and the land sur- 
face will be reduced. Conversely, if the unstable condition 
occurs, the turbulent mixing is intensified and exchanges 
between the air and the land surface will increase. The 
true neutral case is rare. Thus, inclusion of the stabili- 
ty effects provides a better estimation of heat and mois- 
ture fluxes. 

Verification of a large-scale hydrologic model such as 
this is extremely difficult because data are not generally 
available. Some of them are even not well understood. In 
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order to understand the effects of the parameters on the 
GHM, a wide range of sensitivity study is justified. The 
sensitivity studies that were performed are described as 
follows . 

1) Different initial and boundary conditions. 

Natural hydrologic conditions vary temporally and spa- 
tially. This is also true for initial and boundary condi- 
tions in the GHM. For example, the moisture conditions in 
the rainy season are totally different from those in dry 
months. Also, the moisture content in a high rainfall re- 
gion is different from that in a desert. This situation 
requires a great deal of data to prescribe the initial and 
boundary conditions everywhere at a given time. Since 
these data sets are not likely to be available in the near 
future f/ the problem of how long and to what extent the ef- 
fect introduced by approximate or arbitrary initial and 
boundary conditions will persist should be studied. Two ex- 
treme moisture initial conditions (wet and dry) and two 
currently used lower boundary conditions are studied. 

2) Selected physical parameters. There exists a wide la- 
titude of variation in the parameters" that are used in the 
GHM. Many of them can not be prescribed precisely for the 
large-scale parameterization. For example, transpiration is 
influenced by canopy resistance and the distribution of the 
root system, but there are not sufficient data available 
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for the global characterization used in the model. The 
sensitivity study can show us to what extent the effects 
from the variation of each parameter are and then can pro- 
vide some insight for improving GHMs in the future. 

The ground hydrologic model is designed to interact with 
a GCM, but for the purpose of model development and testing 
it is much more practical and efficient to conduct the 
study non-interactively . This means the GHM is driven di- 
rectly by prescribed atmospheric variables . in the PBL that 
are derived from the GCM output data. For this study, the 
input data is available globally for 45 days from July 10 
to August 25 in 1975 from the NASA GLAS GCM. Only eight 
regions with different climate and different land surface 
characteristics across the North America continent are se- 
lected/ for the model evaluation and the sensitivity stud- 
ies. The eight regions, following the code index (I,J) in 
Alfano (1981) and Tsang and Karn (1975) are listed in Table 
1.3*1 with * mark. The relation between (I,J) and (longi- 
tude, latitude) is shown in Table 1.3*2 where the number in 
the parenthesis on the upper left indicates the longitude 
and the number on the lower right indicates the latitude. 
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Figure 1. 3. 1 Measured mean-annual water and energy 
balances of adjacent grass-covered and forest-covered 
catchments, in central Wales, U. K. , (Mintz, Sellers 
and Willmont, 1983) 

































Chapter II 

SPECIFICATION OF LAND SURFACE CHARACTERISTICS 

Land surface hydrological processes are strongly influ- 
enced by the properties of land surface cover and the soil 
underneath. Although land represents only around 29.2# of 
the earth's surface (146*10^ km^) , its effects on the at- 
mospheric general circulation and climate become signifi- 
cant when air moves over the continents. The land surface 
cover is never uniform. It varies spatially from complete 
bare soil to fully covered vegetation. It varies dynamical- 
ly according to seasonal variation of vegetal surface cover 
due to .climatic adjustment and biological evolution of 
plants.'- So, an important, prerequisite step to the devel- 
opment of the GHM is the characterization of soil and vege- 
tation from grid to grid and with seasonal changes. 

Despite the fact that a one-dimensional transport of mo- 
mentum, heat and moisture is applied to each grid, the hor- 
izontal variability of the GHM is accomplished by varying 
the land surface character from grid to grid. In this 
chapter, this necessary information about the global land 
surface will be provided, much of which has been described 
in Lin et al. (1978) and Alfano (1981). In order to real- 
istically model land surface processes spatially and tempo- 
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rally, the following characterization should be provided to 
every grid in the GHM: 

1) specification of soil types and their relevant proper- 
ties. 

2) specification of vegetation types and their relevant 
properties 

3) specification of vegetation fraction and its seasonal 
variation 

2 . 1 SOIL TYPE AND ITS DISTRIBUTION 

In this subsection, the available global soil data will 
be evaluated in order to supply the GHM with appropriate 
information . 

Soil consists of solid phase (organic matter and miner- 
al), liquid phase (water) and gas phase (air, water vapor 

\ 

and other gases). The size of mineral particles ranges 
from very fine clay particles of less than 0.002 mm in di- 
ameter to coarse sandy particles of up to 2 mm in diameter. 
The U.S. Department of Agriculture, based on the fraction 
of sand, silt and clay contained in the mineral matter of 
soil, has presented a classification of 12 soil textures 
(Bridges, 1978) which is shown in a textural triangle in 
Figure 2.1.1. The classification is critical because the 
water movement inside the soil is dependent upon the tex- 
tural composition. In the GHM, 


following Alfano's work, 
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the soil is characterized by its texture into 5 types that 
are consistent with the 12 soil textures in Figure 2.1.1. 
The five types of soil labeled as A to E are listed below : 
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Soil type 
label with letter 

Class by 

textural specification 

A 

Sand 

B 

Sandy loam 

C 

Loamy sand 

D 

Loam 

E 

Silt clay loam 


The percentage of sand, silt and clay in the mineral for 
the 5 types of soil was obtained by Buckman and Brady 

(I960) and is shown in Table 2.1.1, 

/ 


Soil type 

% of sand and 
silt in mineral 

# of clay 
in mineral 

Sand (A) 

95 

5 

Sandy loam (B) 

90 

10 

Loamy sand (C) 

95 

5 

Loam (D) 

65 

15 

Sil.t clay loam (E) 

TO 

3° 


Table 2.1.1 Percentage of soil, silt and clay 
in the GHM soil type 
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Soil Types in GU M 
A-l Sand, A-2 Silt 
B Sandy loam 
C Loamy sand 
D-l Loam, D-2 
E— 1 Clay loam, 

E-3 Silty clay loam 


Soil Types Not in ohm 
S- l Clay 
S-2 Sandy Clay 
S-3 Silty Clay 


Silt loam 

E-2 Sandy clay loam 



Figure 2.1.1 


1 textural soil representation of 
the U.S. Department of Agriculture 
textural class triangle. 
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The water content in tne soil is an important variable 
for determining matrix potential, hydraulic and thermal 
properties. Different types of soil have different water 
holding capacity. 

Ritjema (1970) grouped the soil into 20 types and pre- 
sented the relationships between soil moisture and the soil 
properties for each type-of soil. Ritjema’s classification 
was related to the classification with 5 soil types (Alfa- 
no, 1981) for use in the GHM. Table 2.1.2 shows the rela- 
tion and gives the water holding properties for different 
types of soil. In the table (Alfano, 1981), 0^ c is defined 
as the value when the matrix potential =-350 H^O-cm and 
0p W p is defined as the value when ^=-1 6 ,000 H 2 0-cm. 

The five types of soil in the GHM should be related to 

the soil distribution (Eyre, 1968) in the world to specify 

\ 

the soil type in every grid. Figure 2.1.2 (Strahler, 1963) 
shows the map of the world soil group distribution. The 
linkage between the soil type in the GHM and the zonal soil 
categories (Mitchell, 1976; Chow, 1964) which then can be 
related to the world soil groups (Strahler, 1971) was 
schemed by Alfano (1981) and is used in this study. (see 
Table 2.1.3) 

The previously defined relationship between zonal soil 
categories and the 5 textural soil types permits a global 
distribution of soil to be specified on the basis of the 
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zonal soil categories. Alfano has presented the distribu- 
tion in a map form as given in Fig. 2.1-3a and b. 


/ 
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Soil Type 

vO- 

fc 

0 

pwp 

6 

avm 


Coarse Sand 

.020 

.012 


• 395 

A Medium Pine Sand 

.066 

.023 

.043 

• 35 

Find Sand 

.14 | 

j 

.04 

i 

.10 

• 36 

Humors Loamy Medium 
Coarse Sand 

H 

. 1 1 

.22 

.47 

C Light Loamy Medium Coarse 
sand 

.22 

.10 

.12 

.40 

Loamy medium Coarse sand 

.16 

.02 

.14 

• 30 

Loamy fine sand 

i 

• 13 

.06 

.07 

. 44 

1 

B Sandy loam 

.16 

.06 

.10 

.47 

Loose loam 

.26 

. 1 1 

.15 

. 46 

Pine sandy loam 

.25 

.09 

.16 

• 50 

Silt loam 

• 31 

.09 

.22 

• 51 

D Loam 

.27 

.10 

.17 

.50 

Sandy clay loam 

• 30 

.18 

.12 

.43 

E Silty clay loam 

• 32 

.19 

.13 

.48 

Clay loam 

• 38 

.26 

.12 

• 45 

Light clay 

• 33 

.21 

. 1 1 

.45 

Silty clay 

.44 

.26 

.18 

• 51 

Basin clay 

.48 

• 32 

.16 

• 54 

peat 

1 

.27 

.41 

.86 

_ 


Note: The first one in every type is used in the GHM 

Table 2.1.2 Characteristic soil moisture values, 0„ 
a f C 

y pvp and 9g from Rijtema's (1970) soil types. 
















































































World Groups 
( St rah} er ,1971 ) 

Zonal Category (NO) 
( Michell , 1 976) 

Textural soil 
Type (letter) 

podzolized soils 

Podzol (2) 

Gray brown (3) 
Mtn Valleys (9) 

Sandy loam (B) 
Loamy sand (C) 
Sandy loam (B) 

Lateritic soils 

Lateritic (5) 

Silty clay 
Loam (£) 

Grassland soils 

Prairie (4) 
Chestnut brown (7) 
Chernozems (6) 

Loamy sand (C) 
Loam (D) 
Loam (D) 

Soils of Arid 
region 

Sierozgems (8) 

Sand (A) 

Spils of .cold 
region 

Tund ra ( 1 ) 

Sandy (B) 


Table 2.1.3 Linkage among soil type, zonal soil 
categories and world soil groups (Alfano, 1981) 
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Figrure 2.1.2 World map of principal zonal soil aroups (Strahler, 1963) 
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Figure 2.1-3a World zonal soil categories (Kellogg, 1941; Mitchell, 1976), 

western half of the world. 
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2.2 VEGETATION TYPE AND ITS DISTRIBUTION . 

Attention should be paid to the vegetation cover when 
studying land surface processes. A canopy acting as an in- 
termediate interface between the air and soil is more ac- 
tive than a bare soil surface not only because its struc- 
ture can stimulate turbulence to intensify exchange between 
the air and land surface, but also because its root system 
can take up more water from the deep soil layer to maintain 
a high evapotranspiration rate. Also, vegetation can adjust 
the atmospheric demand by its own biological reaction such 
as opening and closing stomatal pores to control transpira- 
tion. The behavior of canopy makes modeling the vegetation 
functions more important and also more difficult. 

The spatial vegetation distribution strongly depends on 

/ 

the climate which forms vegetation zones called zonobiom- 

\ 

ass. This zonobiomass are more or less modified by the soil 
properties, orography and other factors. The classifica- 
tion with eleven types for the world vegetation was used by 
Alfano (1981), based on Strahler's paper. They are: 

1 ) Equatorial and tropical rainforest 

2) - Temperate rainforest 

3) Evergreen hardwood rainforest 

4) Raingreen forest, wood land, scrab and savanna 

5) Steppe and prairie grasslands 

6) Dry desert and semidesert 

7) Summergreen deciduous forest 
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8) Needleleaf forest 

9) Arctic tundra 

10) Icecaps 

1 1 ) Highland areas 

The detailed description of the 1 1 types of vegetation 
can be found in Alfano (1981). The distribution of vege- 
tation was described by' Eyre (1968) and Strahler (1971). 
The distribution of the different kinds of vegetation was 
derived from the best available detailed map for parameter- 
ization of the global land surface, and is typical of other 
distributions (Collinson, 1977)* Figure 2.2-1 shows the 
distribution of vegetation for all the land grids in the 
world (Alfano, 1980) 

The volumetric organic matter fraction in a soil for 

/ 

vegetation types is given in the table below. 



Type 

(Strahler , 1 971 ) 


Volumetric organic 
matter fraction 
(Buckman, I960) 


Rainforest 

15. 

Temperate rainforest 

10. ; 

Evergreen hardwood forest 

7.5 

Raingreen forest 

10. 

Steppe and prairie grasslands 

5.0 

Dry desert and semidesert 

2.5 

Summergreen deciduous forest 

10. 

Needleleaf forest 

7.5 

Arctic tundra 

5-0 

Highland areas 

5-0 
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Figure 2.2-la World vegetation types (Strahler, 1971), western half of the world. 
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Figure 2.2-lb World vegetation types (Strahler, 1971) ; eastern half of the world. 



2.3 SUBGRID VEGETATION PARAMETERIZATION — DISTRIBUTION 

ANTTVA'HlITION ~ 

The vegetation density is represented by the vegetation 
fraction 6 C in a grid and differs from grid to grid and 
varies with season. The definition of vegetation fraction 
is the fractional portion of vegetated cover in a grid. 
Then, the bare soil fraction is 6 ^ = 1-6. So, if 6^=0. 6, 
it means that a grid is covered with 60# of vegetation and 
40# of bare soil. Based on the pictorial difference in 
color and tone from the map (Strahler, 1971; Peele, 1975; 
Lee , 1977), Alfano (1901) assumed that the density corre- 
sponds to the full growing season and derived the global 
distribution map by the minimum bare soil fraction <5^ as 
shown in Figures 2 . 3 - la and 2-3-lb. 

Figure 2.3-1 specifies the minimum bare soil fraction 
distribution for all land grids. However, the vegetation 
fraction changes with seasons. Actual vegetation fraction 
should be modified on the basis of the minimum ^ and the 
annual percentage of the period in foliage. 

To consider the seasonal change of 6 (Figure 2.3.2), 

V 

Alfano makes the following suggestions: 

l) The period of vegetation in foliage, in days, which 
is correlated to the seasonal change in solar energy, con- 
sists of a 17.5 # (Rosenberg, 1974, Money, 1976) portion 
for growth (spring), a 65 # for full foliage (summer) and a 
17.5 # for senescent (fall) shown in Figure 2.3-2 
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?) The foliage development is always centered on the 
summer with the following key Julian dates: 

a) the start of vegetation development (JGSS) 

b) the beginning of the full summer foliage (JGSSP) 

c) the end of vegetation development (JGSPP) 

d) the end of the foliage period returning to winter 
conditions (JGSP) 

e) the number of days in the foliage or green season 
( JDGS) 

3) The annual percentage of the period in foliage 
(JDGS/365) is specified as a linear function of latitude 
determined from the typical foliage pattern (Money, 1976) 
with continual full foliage below 22° latitude and no fo- 
liage change above 66.5° latitude (Lee, 1980) as shown in 
Figure 2.3*3* 
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world (Strahler, 1971; Peele, 1975; Lee 1977). 
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Figure 2.3.2 An example of seasonal variation for New Jersey. 
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Figure 2.3.3 Percentage of vegetation foliage for the year 
as a function of latitude (Alfano, 1981) 
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2.4 SOIL PROPERTIES 

t 

Two kinds of soil properties, the hydraulic property and 
the thermal property, must be considered in the GHM. These 
properties are related to the soil moisture content. 

When one deals with the water movement in soil, the ba- 
sic principle for moisture flux q in the soil is Richards’ 
law (Richards, 1931) which is expressed by (for one-dimen- 
sional) 

q=-K* a '//az + K (2. 4-1. a) 

or =-D*ae/az +k (2.4.i.b) 

D= K*4f/<\Q (2.4.1 -c) 

where y is the soil matrix potential, K is the hydraulic 
conductivity and D is the diffusivity. 

It should .be pointed out first that, from a physical 

/ 

view point, the Eq. (2.4. 1b, c) doesn't add any information 
to Eq. (2.4.1a). On the other hand, the expression in Eq. 
(2.4.1b) will be difficult to treat at the interface of 
heterogeneous soil because there is normally a sharp dis- 
continuity in the soil moisture contents across the inter- 
face of layered soil and only the potential remains contin- 
uous across the interface (Hanks and Ashcroft, 1980). 
However, the expression is more acceptable in the general 
case because of its explicity and simplicity. 
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2.4*1 Total water potential and matrix potential 

The total water potential ti is formally defined as the 
amount of work that a unit quantity of water in an equilib- 
rium soil water system requires when it moves to a pool of 
water in the reference state at the same temperature (Hanks 
and Ashcroft, 1980). Disregarding the osmotic potential 
and pressure potential, consists of matrix potential and 
gravitational potential: 

H ='V / + z 

The gravitational potential is independent of soil na- 
ture and equal to the vertical distance z between a refer- 
ence point and the point in question. 

The matrix potential 'f is related to the adsorptive 

force of the soil matrix (hence, "matrix potential") and 

/ 

equal to the vertical distance between that point in the 
soil and the water surface of a manometer filled with wa- 
ter and connected to the soil point in question via a ce- 
ramic cup. It is a function of soil moisture, soil type and 
so on. 

For a given soil, the matrix potential of moisture in 
the soil is not a unique function of the soil moistfure con- 
tent. This phenomenon is known as hysteresis (Hanks nd Ash- 
croft, 1980). In Figure 2.4*1, the wetting and drying pro- 
cesses produce an envelope that gives the extreme ranges of 
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possible water content that can be associated with any par- 
ticular matrix potential. 

The moisture matrix potential curves used in this GHM 
for five types of soil has been given by Ritjema (1970), 
which are obtained from measurements taken during the soil 
drying process. They are shown in Figure 2.4*2 (Alfano, 
1981 ) • 


/ 




Q [vmc]’ 


Figure 2.4.1 Soil suction curves for desorption and 

sorption as a function of moisture content 
(Hillel, 1977; Slayter, 1977). 
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2.4.2 . Hydraulic conductivity 

The hydraulic conductivity is related to the moisture 
content. The hydraulic conductivity is very difficult to 
determine in the field and can only be measured (Nielsen et 
al . , 1964; 197?) in the field under very favorable condi- 

tions. Even in the laboratory, the experimental measure- 
ment is also difficult (Nilsen and Biggar, 1961). The em- 
pirical formulas suggested by Ritjema (1970) are 

K( 0) = a r *(-f(e )) -1 * 4 •+<% (2.4.2a) 

K(0)= K r *exp(oi r *'V') (2.4.2b) 

where a r , o* r , K r , *h r was given by Ritjema (1970) and shown 
in Table (2.4.1) 

Alfano, based on Eq. (2.4.2), has plotted the curves of 

K(0) versus volumetric moisture content for five types of 

/ 

soil which are shown in Figure 2.4-?. 
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2 . 4 • 5 Diffusivity 

% 

Alfano calculated diffusivity with Eq. (2.4.1c) using 
the values of K and ^ given by Ritjema (1970) which are 
plotted in Figure 2.4-4 for five types of soil. 

In order to make the integration of K and D more conven- 
ient, the relation of K and D with the soil moisture con- 
tent was approximated in several segments by the following 
equations (Alfano, 1981): 

K(©)= a^ exp(a 2 (©-a^) ( 2 . 4 . 4 a) 

D(0)= b 1 exp(b 2 (©-'b^) (2.4*4^) 

where a^ , a 2 , a^ , b^ , b 2 , and b^ are constants which will 

only change the magnitude from one segment to another. The 
constants are listed in Table 2.4-2. 


/ 
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2.4*4 Thermal property — heat conductivity and heat 
capacity 

Heat exchange in a soil is primarily caused by conduc- 
tion. The heat flux in the soil can be described by Fouri- 
er's law. For vertical heat flux in a one-dimensional col- 
umn, it can be expressed by 

G = - AdT/dz (2.4-5) 

where G is heat flux positive downward and A. is conductivi- 
ty (cal cm -1 s - ^ deg"" 1 ) • The conductivity A is expressed as a 
linear function of the volumetric fraction of the compo- 
nents such as sand and silt, clay and organic matter and 
given by 

3 

A= ZA 5 v. (2.4.6) 

i = l 

where is the volumetric fraction for sand and silt, 

clay and organic matter, and is the corresponding con- 

/ 

ductivity. In the GHM, The heat conductivity A^ was ob- 
tained by fitting linear segments to the curves given by 
Sellers (1961) as follows while considering moisture ef- 


fect (Alfano, 1981): 

For sand and silt, 

^ =0.001 (61 0+0. 5) 0<O .05 ( 2 . 4 . 7 a) 

=0.001(14.20+2.85) 0-05<G<0.18 (2.4. 7b) 

=0.001(1.860+5.1) 0> 0 . 1 8 (2.4.7c) 

For clay, 

>- 2 = 0 . 001 (210+0.4) 6<0.14 ( 2 . 4 . 7 a) 

=0.001(5.130+5.0) 0>O • 1 4 ( 2 .4 • 7e ) 
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For organic matter, 

A-3=0. 001 (0.426+0.1 ) 0<O.24 (2.4.7f) 

=0.001(1.70-0.17) 0>O.24 (2.4.7g) 

Another important property for soil is the heat capacity 
of the soil. Sellers (1965) suggested that the soil volu- 
metric heat capacity C can be estimated by averaging the 
heat capacities for the components in the soil: 

C= 0.46(V 1+ V 2 )+0.6V ? +0 i (cal/cm 5 /°K) (2.4-8) 
where (V^l^) is the volumetric fraction of mineral, is 
the volumetric fraction of organic matter and 0^ is the 
volumetric moisture content in layer i. 

2-5 VEGETATION MORPHOLOGY AND PHYSIOLOGY 

With regard to heat and moisture transfer in a soil- 

/ 

plant-atmosphere continuum, the canopy morphology and phys- 
iology play a very important role. For the canopy morpholo- 
gy, there are four important parameters: canopy height h , 

V 

zero displacement thickness d and roughness height z Q re- 
lating to the turbulence flow in the atmosphere, leaf area 
index LAI relating to light penetration from the canopy top 
down to the vegetated ground surface, and the root system 
in soil which relates to the water uptake. 

For the canopy physiology, the most important factor for 
hydrologic modeling is the canopy resistance which adjusts 
the transpiration rate from leaves. 



57 

2.5*1 Vegetation morphology 

t 

In the GHM, the whole canopy is treated as a bulk layer. 
Only the integral effect of the localized parameters is 
considered in the bulk layer model. 

1) Leaf area index LAI: 

Leaf area index is defined as the total area of leaves 
over unit area of the ground surface. Its magnitude influ- 
ences heat and water transport processes. 

First, when the light penetrates from the canopy top 
down through the canopy, most of it is intercepted and ab- 
sorbed by the leaves. The residual, which can arrive at the 
ground surface, may be approximately expressed by an expo- 
nential decay law: 

< R n (bot) = R n (top)exp(-X c *LAI) (2-5-1 ) 

\ 

where R n (top) is the net radiation at the top of the cano- 
py, R n (bot) is the portion of the net radiation reaching 
the ground surface, X c is the extinction coefficient which 
depends on the type of canopy, the leaf surface orientation 
and the solar elevation angle. LAI is the leaf area index. 

Second, the absorbed net radiation by the canopy will be 
partitioned into latent heat, sensible heat and heat stor- 
age inside the bulk canopy layer. The heat storage capaci- 
ty depends on the heat capacity of the bulk layer which is 
primarily determined by biomass and moisture contained in 
the canopy. 
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Third, the actual transpi ration • is dependent on the can- 
opy resistance and is expressed by 

E c = ^a*^* {T c ) - < la^( r v +r cr ) (2. 5*2) 

where q (T ) is the saturated specific humidity, q_ the 

specific humidity in air, r v the air resistance to the ac- 
tual transpiration and r cr the canopy resistance to the ac- 
tual transpiration. The- canopy resistance is approximately 
related to the leaf resistance by (Landsberg, 1975) 

r c = r-j/LAI (2.5-3) 

where r^ is the leaf resistance for the canopy which is 
approximately equal to the stomatal resistance r . 

There are many articles reporting leaf area index val- 
ues for different kinds of plants (Landsberg and Cutting, 
1975; Montheis, 1975; Perrier , 1 981 ) . Prom the available 
data it; was found that the value of leaf area index depends 
on the species. Table 2-5 -1 exemplifies the values of LAI 
for different plants (Perrier, 1981). In the GHM, an aver- 
age value LAI = 3 is taken for a short canopy such as 
grassland and crop, and LAI = 6 or 7 for a tall forest. 
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Plant 

Leaf area index (LAI) 

Evergreen forest 

7-10 

Bereal coniferous zone 

2-10 

Mediteranean semipervirent forest 

4-5 

— 
Mozon forest and mixed dry forest 

6-10 (even to 15) 

Decidous temperate forest 
(total vegetation) 

CO 
vo 1 

1 VO 
Ttw 

Woody savanna 

1- 1.5 (dense) 

.2-. 4 (dry weather) 

2- 4 (rain season) 

Temperate grassland 
Douglas fir 

5-6 

-6 

(projected area base) 
(Landsberg, P71 , 75) 

Grassland 

4.4 

(Landsberg, P44, 75) 


Table 2.5.1 


Leaf area index for different canopies 
( Perrier , 1 981 ) 
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2) Root density distribution 

The transpiration from leaves to the air depends in part 
on the ability of the root system to absorb water from the 
soil. The uptake capacity is dependent on the root distri- 
bution and development. Available data for the root distri- 
bution in soil is not sufficient to give a systematic sum- 
mary or quantitative description of root. Its field 
measurement is very difficult. There are some scattered 
data for the crop root system but only a few are related to 
the forest root system. (Miller, 1938; Ritjema, 1970; Hut- 
tel, 1975; Hillel, 1977; Wood, 1980; Perrier, 1981; Waggon- 
er and Turner, 1 97 1 ) - 

The root density should be properly defined as the den- 
sity of root surface because the root weight may not be 
correlated with the mechanism of the water uptake in for- 
est. The root density in the first layer is defined as 
equal to the surface area of roots in layer 1 divided by 
total surface area of roots in the whole root zone. 

Table 2.5.2 lists some of the available data for evalua- 
tion in the GHM. Figure 2.5-1 shows the distributions of 
some root systems. . 

In the GHM, the distribution used is the same as the 
shallow and deep root systems suggested by Hillel (1977). 



Canopy 

Root density 
R at 

in top 10 cm 

Root density 
from 10 cm 

to around 1 m 

Reference 

Ivory 

forest 

0.21- 0.47 

• 0.79-0.53 

Huttel (75) 

Hardwood 

woody 

0.50 

0.50 

Wood (80) 

Savanna 

0. 50 

0.50 

Perrier(81 ) 

Grass 

land 

/ 

0.50 

0.50 

Perrier(81 ) 

Shallow root 
system 

0.50 

0.50 

Hillel (77) 

Deep root 
system 

0.25 

0.75 

Hillel (77) 


Table 2.5*2 Root density distribution 
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2.5*2 Canopy physiology — canopy resistance 

Canopy resistance r , in general, may be considered as 

the impedance to water transport from the vegetal surface 

to the air. r c is approximately equal to leaf resistance r^ 

divided by leaf area index (LAI). The leaf resistance r-^ is 

almost equal to the stomatal resistance r provided the 

s 

stomata are not closed entirely (Monteith, 1975) 

Canopy resistance (or leaf resistance or stomatal resis- 
tance) is affected by many factors such as light, wind ve- 
locity, air temperature, air humidity and leaf water poten- 
tial (Landsberg and Cutting, 1975) (see Figure 2.5-2). 

In this study, Eq. (2.5*2) is not used for calculating 
actual transpiration because of uncertainty of the canopy 
resistance r . The actual transpiration is calculated by 

y v r- 

scaling from the potential (or % maximum) unstressed transpi- 
ration E u which is defined as follows: 

E u ■ VC<l*(' l; o>-'laJ''( r v +r o n ’ <2-5-4) 

where r Cm is the minimum unstressed canopy resistance. In 
all succeeding sections, r is used to denote this minimum 

V 

resistance r . The minimum canopy resistance is assumed 
cm . 

to be a constant which depends on the vegetation type. r c 
= 75 s/m for tall canopy and r c = 100 s/m for short canopy 
are used in this study. 
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ORIGINAL FAGS IS 
OF POOR QUALITY 



Figure 2. 5-2b Relationship between stomatal 
resistance r s and vapor pressure deficit 6 e for 
three ranges of soil matrix potential r tl 
for Douglas Fir (Landsberg and Cutting, 1975) 
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Figure 2.5-2c Wheat crop resistance versus net radiation 
(a), air deficit (b) and wind speed plotted 
as a function of the crop aerodynamic 
resistance (c) (Perrier, 1981). 
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2.6 ALBEDO 

Albedo is defined as the fraction of shortwave radiation 

reflected by the land surface cover. According to this def- 

inition, the net radiation absorbed by a surface is 

R n = R s*< 1 - a l> + R LN (2.6.1) 

where R g the shortwave radiation 

R LN — t ^ ie net lon g wave radiation (amount absorbed 
minus the amount emitted by a surface). 

a^ albedo (in percentage) 

Albedo depends on several factors such as sun elevation, 
direct solar to scattered radiation ratio, land surface 
cover and surface conditions (moisture, roughness and so 
on). Kondratyev (1981) and Carson (1981) have discussed the 
available information in detail and presented the parame- 
terization of albedo for use in numerical climate models. 

In the GHM, Alfano's scheme is followed, which was based 
on the work done by Idso et al. (1975)» Strahler (1971), 
Monteith (1975) and Pederer (1968). The albedo of soil is 
a function of soil moisture. It can be expressed as a seg- 
mental linear function and given by 

a lb =30, " 2 -* e avf /0 ' 6 0<9 avf <0.6 

=28.-l6.*(e avf -0.6)/0.4 o.6<e avf <i • 

where © avf = “ e pwp)/( 0 fc" & pwp) and e l is the volu “e tr i c 
moisture content in the surface layer. 
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The canopy albedo a^ c is determined by the canopy type 
and the water stress condition (or soil moisture) as well 
as the location where the canopy grows (Posey and Clapp, 
1964)* Corby et al. (1977) gave the albedo distribution 
along the latitude for use in a 5-level atmospheric model. 
In the GHM, the albedo for the different canopies is shown 
in Table 2.5*3 (Alfano, .1981), which accounts for canopy 
type and soil moisture content. 


/ 
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NO. 

Vegetation type 

Albedo 

1 ,2 

Tropical & Temperate 
Rainforest 

14 

3 

Evergreen hardwood 
Forest 

12 

4 

Raingreen forest 

20 - 8 *( 0 -e pwp )/(e fc -e pwp ) 

5 

Grasslannd 

20 

6 

Desert 

30 - 1 5* ( S-» pwp ) / ( 0f c - e pw p ) 

7 

Deciduous forest 

18 

8 

Needleaf forest 

18 

9 

Arctic tundra 
1 *JDAY 125 
1 26 JDAY 125 

Spring ( NH) : 39-1 9( JDAY/ 1 25) 
Fall (NH) :20-19(JDAY/l25) 

\ 

10 

Ice 

40-80 (Tsang & Karn, 1972) 

11 

Highland Areas 

25 -io*(0-9 p „ p )/(e fc -e pwp ) 

12 

OCEAN 

7-9 (Tsang & Karn, 1972) 


0 = O.20 i + 0.80,, 

JDAY: Julian day 

NH: Northern Hemisphere 


Table 2.5*3 Vegetation albedo parameterization 




Chapter III 


DESCRIPTION OP MOISTURE MOVEMENT AND HEAT 

TRANSFER 

Only the one-dimensional vertical transports of the 
moisture and heat are considered. The exchanges of water 
and heat are between the air and soil, the air and canopy, 
the canopy and vegetated ground surface, and the soil in 
r the root zone. Figures 3*1-1 a,b show the schematic of the 

processes. In the figures, the subgrid parameterization is 
accomplished by dividing a grid into two parts : canopy and 
bare soil. The bare soil fraction is <5^ and the rest of the 
grid, 6 c = (1.- ^), ia the canopy fraction. For moisture 

movement, the root zone is divided into two layers : a thin 

/ 

surface layer over a lower deep layer extending to the 
bottom of the root zone. For heat transport, a surface 
layer has a thickness of ' 1 cni and the whole layer thickness 
(see Figure 3*1 -1b) equals the depth of penetration of the 
annual temperature wave* The thicknesses of the layers 
used for the soil moisture movement do not coincide with 
those used for the heat transport. Over the vegetated 
ground surface, one more canopy layer is added. 


- 71 
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Figrure 3.1- la 


Moisture movement 




Ficnire 3.1-lb Heat transport 
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The time-dependent variables are the soil moisture con- 
tent 0^ in the surface layer, &2 in the lower layer, the 
bare soil surface temperature T the vegetated ground 
surface temperature T gc , the bulk canopy temperature T c , 
and the layer average temperature T* of the whole layer. Ac- 
companying outputs are the grid average evaporation and 
sensible heat flux H^. (note: later all variables preceded 
by 'grid average' are called 'grid'). 

In order to drive the GHM, the atmospheric variables are 
the temperature T , the wind velocity u m , the specific hu- 
midity of air (subscript m means the mean value in the 
PBL) , precipitation P r , air pressure P and net radiation 
R n - They are generated from the GCM (Tsang and Karn, 1973; 
Halem etal., 1979) results. 

The ''GHM deals with a dynamic system which, from a mathe- 
matical view point, is an initial-boundary value . problem. 
The upper boundary conditions are obtained from coupling 
with the PBL, which determines evaporation and sensible 
heat flux in the surface layer of the PBL. The initial and 
lower boundary conditions should be prescribed. The formu- 
lation for the GHM is presented belbw section by section 
based on the physical phenomena. 
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5-1 GOVERNING EQUATION FOR SOIL MOISTURE 

The saturated volumetric moisture content 0 is defined 

8 

as the volume of water per unit volume of soil when the 
soil is fully saturated. 

The moisture content at field capacity 0 fc is the volu- 
metric moisture content that a soil reaches and maintains 
after it has been thoroughly wetted and allowed to drain 
freely for a day or two. 

The moisture content at permanent wilting point 0^ is 
the volumetric moisture content when water extraction by 
plants has almost ceased. 


The definitions of the moisture content at field capaci- 
ty and permanent wilting point are idealized and difficult 

/ 

to determine uniquely for different soils and different 
plants. In order to overcome this difficulty, the defini- 
tions used in this study are the quantitative ones given in 
Chapter 2. e fc is defined as 0 at which the matrix poten- 
tial equals to -350 H 5 0-cm. 0„ 1JT% is defined as 0 at which 

L. PWP 

the matrix potential equals to -16,000 HgO-cm. 


Maximum available moisture content 0 Is 

avm 

®avm = ®fc” ®pwp 

Available moisture content 0 is 

av 


®av ~ ®pwp 

The fraction of available moisture content 0 av ^ is 


®avf = ®av^®avm 



The soil moisture movement in the two layers are sche- 
matically shown in Figure 3. 1-1 a. Without regard to the 
rain interception by canopy and intercell exchange, the 
variation of volumetric moisture content 0^ in the surface 
layer for the GHM is controlled by the following compo- 
nents: the surface infiltration 1^ (into the bare soil sur- 
face) and I c (into the vegetated ground surface), the eva- 
poration from the bare soil, the water uptake by the 
roots in layer 1, the interfacial moisture flux q^ ^ between 
layers t and 2 and the subsurface runoff when moisture 
content 0^ exceeds the field capacity. Storage change of 
liquid water above the ground surface is neglected in the 
GHM. The equation for 0^ is expressed by 

d, de, /dt=( V 0 ! ) V'll 2- R S 1 +I 81 

. V' • ■ ■■-, >•. -• •• • - (3.-1.1) 

where I a ^ is the infiltration from melting snow and ice and 
is disregarded in this study. 

The volumetric moisture content for the second layer 0 2 
is governed by 

^2^®2^ ^ = ” 2 — ^ , 23""^s2"*" ^s2 (3*1 *2) 

where U2 is the water uptake by the roots in layer 2, q2^ 

the interfacial moisture flux, R„„ the subsurface runoff 
from layer 2 when the moisture content in the layer exceeds 
the field capacity and 1^ is the contribution from ice 
melting and is also disregarded in this study. 
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The moisture content in the third layer below the root 
zone is not known. But it seems likely that the soil mois- 
ture condition in the third layer changes very slowly with 
time and may be expressed by 

. 0^= constant + f(£t) (£<<1) (3.1*3) 

However, for the short period experiments, 0^ may be ap- 
proximated by 0£ c for this study. Another lower boundary 
condition is also used as part of the sensitivity study. It 
is a constant flux at zero moisture gradient condition 

^23 = ^23 5©/ez| • (3*1 *4-) 

where K 2 j is defined in Eq. (3*3*4-) 

In Eqs. (3-1.1) and (3*1*2), the parameters d^ and 
are determined based on the root zone thickness and root 
system distribution. For most plants, the top 10 cm below 

the surface has a denser root distribution' and is selected 

/ 

as the surface layer. The deep root zone thickness depends 
on the vegetation type. In this study, a shallow root sys- 
tem and a deep root system (Hillel, 1979) have been used 
for the short canopy and the tall canopy, respectively. 
They are d^ = 40 cm for the short canopy and d 2 = 90 cm for 
the tall canopy. Different values of d 2 are used in the 
sensitivity study. 

The fluxes in Eqs. (3*1*1) and (3*1*2) are treated in 


the subsequent sections. 
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3*2 GOVERNING EQUATION FOR GROUND TEMPERATURE 

Figure 3«1-1b shows the scheme for the heat transport 
between the air and a bare soil surface, the air and a 
bulk canopy layer, a bulk canopy layer and a vegetated 
ground surface, and in the soil. Sellers (1965) pointed out 
that the soil temperature has layered characteristics: a 
thin surface layer influenced by the diurnal temperature 
cycle and a deeper layer (the annual layer in Figure 
3*1 -1b) influenced by the annual temperature cycle. In the 
surface layer, the ground temperature has a high tempera- 
ture gradient (Bruce et al. , 1977; Kimball et al., 1976). 

Bhumralkar (1975) and Blackadar (1976) have, based on 
the simple harmonic solution of the temperature wave, inde- 
pendently proposed a force-restore method to predict the 
soil temperature variation. According to this method, the 
ground surface temperature is dependent on the heat flux on 
the surface as well as on the restoring influence of the 
deeper soil layer. Deardorff (1978) compared the methods 
with those currently being used in the GCM model and showed 
that the force-restore method can provide a better pre- 
diction. Lin (1980) revised the surface temperature formu- 
lation, which assumed surface temperature to be a linearly 
average temperature within a thin layer of the order of 1 
cm. The final rate equations for the bare soil surface 
temperature T ^ and the vegetated ground surface tempera- 
ture Tg c in the GHM are expressed as 
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(3.2.1a) 
(3.2.1b) 
(3.2.1c) 

tn Eq. (3*2.1) G^ is the heat flux into the bare soil 
surface, G is the heat flux into the vegetated ground sur- 
face, C is the heat capacity of the soil and d 1 is the 
damping depth of the diurnal temperature wave equal to (T 
'V(ftC)) 1 / 2 and T=86400 seconds 

r The layer average temperature f is governed by a rate 

equation 

dT/dt=( <5^G b + <5 c G c )/(C( 365 Tt ) l/2 d“) (3-2.2) 

1 /2 ~ * 

where the term of (365^) * s the penetrating depth of 

the annual temperature wave (Lin, 1980). 

For <the bulk canopy layer temperature T , the direct use 

c 

of heat conservation principles will give 

C c h c dT c /dt = X M (3.2.3) 

where C c is the volumetric heat capacity of the canopy, h c 

is the height of the canopy and X^j is the increment rate of 
the heat storage in the bulk canopy layer of unit sectional 
area. 

Eqs. (3*2. 1 )— (3*2.3) can be integrated if initial condi- 
tions are given. The determination of G^ , G Q and will be 
discussed in the next section. 


oidT gb /dt = 20^(0^ )-27C/r(T gb -T) 
o«dT gc /dt = 2G c /(Cd‘ 1 )-27T/r(T gc -T) 
6L = ( 1 +25/ d 1 ) , S = 1 cm 
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3.3 FORMULATION OF VARIOUS FLUXES OF VERTICAL MOISTURE 

MOVEMENT 

In the previous section, there are many fluxes left to 
be determined. These are soil moisture fluxes q^ between 
layers i and j , subsurface runoff R g ^, evaporation or eva- 
potranspiration from bare soil or canopy, water uptakes U^ 
and U 2 from the root zone, infiltrations 1^ and I c , heat 
conductions G^ and G c and heat storage increment rate X^. 
The evaporation and evapotranspiration are related to the 
sensible heat fluxes from the bare soil, H c from the 
canopy and H ^ from the vegetated ground surface to the 
canopy layer. This section and the next section will be 
devoted to the formulation of these fluxes. 

3.3*1 Moisture flux in soil 

/ 

According to the Richards law, the moisture flux is 
expressed by 

q^ = -Dd©/dz + K (3-3-1) 

where D and K are functions of moisture content 0. Since 
only average soil moisture content 0^ and © 2 are available 
in the GHM, q^ between layers i and j is approximated by 
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z j 


1/d ij 

J 

q ij d2 = - 

. 

DdG/d A j + 

_ - 


0, 




Kd©/ ( 0^-0^ ) (75-3.2) 


Then, the average q^ is 

q ij = D U (G J " 0 i )/d U + K U 

0, 


(3.3.3) 


where = 1 / ( 0^ ) 


Dd© 


0. 


r e j 


and 


and 


hi = 


KdQ 


(3.3.4) 


e, 


d^ is the transport distance between 0 = 9^ and 9 = 

0., depending on the layer thicknesses and the profile of 
J 

soil moisture distribution in the layers. In this study, 
d„ are approximated by 


0 , ) ^ 0 ^ 


/ Vj r WO ' 

d 4i = 0.5(d i + d.) if 

3 3 0 1 < © 2 < 

©. < 0 P > 

= 0.25(di + dj if 
3 3 9 1 > 0 2 < 


0, 


(3.3.5) 


(3.3.6) 


Since the relation between K, D and 0 has been approximated 
with straight-line segments on semi-log plots by Eq. 
(2.4.4) 

K(G)= a 1 exp(a 2 ( G-a^) ) (2.4-4a) 

D( 0) = b ie xp(b 2 (G-b 3 ) ) (2.4.4b), 

the expression for K,. and D i . will depend on where 0- and 

0. are located. Assuming 9^ < 0. without losing generality, 
J A J 



there are three possible cases shown in Figure 3.3.1 where 


0^ b and 02^ are the values of the moisture content at the 
boundary demarcation points of the segments shown under the 
O-term in Table 2.4.2. 

If case A) is true, (D^, K i j) is equal to 

A i j = (A(e i )-A( 0 j)J/ln[A( 0 i )/A(Q ;J )] (3-3.7) 

If case B) is true and 0^^ and 0^ are located in either 
side of a demarcation point 0^ (i=1,2), then 

A ij =1 ^ e j- 0 i ) {((A(0 ib )-A(0 i ))/ln[A(0 ib )/A(0.)] 

*( e ib“ 0 i) + [ A ( 0 j)~ A ( 0 ib)J/ ln (A(0j)/ A (©ib)J[ 9 j- e ib}j 

(3.3.8) 

If case C) is true, 0.^ and ©j are separated by the de- 
marcation points, thus 

< A i j = f/(0 ;j )-0 i ){C(A(0 1b )-A(.0 i )]/ln(A(0 1b )/A(© i )J* 

(0 1b -0 i )+CA(02 b )- A (ei b )J/ 1 nC A (02b)/ A (0i b ) ]* 

(0 2b- 0 1b ) + ( A(0 3 ) - A(O 2b ) J /ln i A(0 j )/A(e 2b ) 3 ( 0 j" G 2b ) J 

(3-3.9) 
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3*3*2 Infiltration at the ground surface after rainfall 

Infiltration is defined as the entry of water into the 
soil body through the ground surface. It is an important 
hydrologic process because its rate determines the amount 
of water which accounts for surface storage and runoff af- 
ter rainfall. 

Given boundary and initial conditions, the rate of in- 
filtration can be calculated with great detail by numeri- 
cally solving the governing equations (Wang and Lakshmina- 
rayana, 1968; Whisler and Bouwer, 1970). However, the 
numerical solution is rather complex and the required soil 
data are also difficult to obtain. This method is not suit- 
able for the GHM. As discussed in the introduction, several 
empirical methods are commonly used to model infiltration 
(Mein dnd Larson, 1973)* In this study, the Green-Ampt ap- 
proach is adapted for the GHM. According to their classic 
paper, infiltration into an initial unsaturated soil gener- 
ally occurs under the combined influence of gravity poten- 
tial and matrix potential. The original formulation (Green 
and Ampt, 1911) for infiltration rate was derived for water 
with a ponding surface into a deep homogeneous soil column 
with a initially uniform moisture distribution. It is as- 
sumed that water enters the soil creating a saturating flow 
zone with a sharp wetting front which separates a saturated 
zone from an unsaturated zone (Figure 3*3*2). 
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Figure 3*3.2 Infiltration 


Omitting the derivation, the infiltration rate I ( I fe or 
I c ) can be expressed as 

I Q = K a (S f + L f + H 0 )/L f (3-3.10) 

where K s is the soil hydraulic conductivity at saturation, 
H q is the depth of water ponding on the surface, is the 

effective matrix potential difference across the wetting 
front known as suction head and is the distance from the 
surface to the wetting front. 

In the GHM, the water ponded on the surface is ignored 
and thus H Q is equal to zero. The infiltration rate I (1^ 
or I ) is determined by the following formulation: 
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f 


I o' = Min(K 8 (S f +L f )/L f , P r /j> w ) (3.3.11a) 

L^. can be determined by 

dL f /dt=I o /(0 a -0j) (3.3.11b) 

where 9^=©^ if <d^ or ©^=©2 ^ L f >d i * 
or =0 if rain stopped or just begun. ( 3 - 3 . 1 1 c) 

All variables in these equations are predictable except the 
suction head S^. Mein and Larson (1973) suggested that the 
suction head can be expressed by 


S 


f 


1 

r 


'Y' • dK r 


K r (0 


in ■ 


where K r = K/K g . 


(3.3.12) 


If the matrix potential and conductivity are given in 
terms of soil moisture, S^. can be calculated directly. 

/ 

3.3.3 Surface and subsurface runoffs 


Surface runoff and subsurface runoff account for the wa- 
ter loss from the soil as shown in Figure 3 . 1 -la. When rain 
is heavy, this loss is important for the water budget. 

In the GHM, the parameterization used by Alfano (1981) 
is followed. In a GCM, the values of grid-average intensi- 
ty of precipitation is in general smaller than the infil- 
tration, and the ponding and evapotranspiration from the 
ponded water can be ignored. Thus, the surface runoff is 
considered as the residual of the precipitation minus the 
infiltration. The formula for the runoff is 
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/ 


R q = Max-(P r /ji w -^I b -^I c , 0) (3-3.13) 


For the subsurface runoff R gi (1=1,2) a general rela- 
tionship, as noted by Sellers (1965) and suggested by Ara- 
kawa and Mintz (1972), was used in this study. R . can be 

3 1 

expressed by 


R gi Dt 


= 0 

= di(0 f 0 fc )r 

= d i ((e g -0 fc )r m+ (0 i -© 3 )) 


0i<0fc ( i=1 ,2) 


Q fc <0 i <0 s 

0 s <G i 


(3.3.14) 


where Dt is a half-hour and r is a scaling factor. Alfano 
(1981), by using observed data (Bruce et al., 1977), sug- 
gested that r value be calculated by 
r=r m ((0i-0 fc )/(S s -0 fc ))“ 

where r =0.5 and m=1 .0. 
m 

/ 

5.3.4 evapotranspiration 


Evapotranspi ration represents the moisture loss from the 
land surface through evaporation from bare soil and ponded 
water, and transpiration . of plants. It plays an important 
role in linking the water transport and energy, transport 

0 t 

process because latent heat flux is supplied when water 
evaporates. The estimation of evapotranspiration is diffi- 
cult and has been studied extensively by multi-disciplinary 
investigators because the process involves the complex in- 
teraction of the atmosphere, plant and soil. There are many 
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empirical formulas that exist for estimating evapotranspi- 

/ 

ration (Fritschen, 1981), among which the combination meth- 
od or Monteith’s formula (Monteith, 1975) has been widely 
used . 

One of currently-used methods is to calculate the actual 
evaporation or transpiration through scaling potential eva- 
poration or unstressed transpiration. 

In the earlier version of GCMs , the concept of the scal- 
ing factor developed by Budyko (1965) bas been extensively 
used, i. e. 

? ■ E / E p 

In the GLAS GCM, the scaling factor depends on available 
water © avf . defined in Section 3*1 in the root zone of 1 m 
thick. In the Lin at al. GHM, the scaling method was re- 
fined. <In each grid, the total evapotranspiratin is the sum 
of the 3ubgrid portions of the bare soil evaporation and 
transpiration from the canopy. The bare soil scaling factor 
^ and canopy scaling factor used to estimate these actu- 
al fluxes from potential evaporation and the unstressed 
tanspiration are adapted in this study. Potential evapora- 
tion indicates the evaporation from A fully saturated soil 
surface with actual ground temperature. For the direct 
coupling with the PBL, the potential evaporation from bare 
soil can be expressed as 

E p = C T u»P a [,*(I g )-q m )] 


( 30 . 15 ) 


89 

where C^, is the heat transfer coefficient derived in Sec- 

* . 

tion 3.4*2, q (T ) is the saturated air specific humidity, 

o 

q ffl is the air specific humidity in the PBL and u* is fric- 
tion velocity. C^u # may also be considered as the inverse 
of -aerodynamic resistance r y to heat transport between the 
ground surface and the location where q ffl is defined in the 
PBL. The unstressed transpiration implies the transpira- 
tion from leaves when plants are under no water stress con- 
dition. The unstressed transpiration may be expressed in 
terms of a canopy resistance as 

E u = J > a(l*( T o ) - , ln.J/( r v +r ou 1 ) (50-16) 

where q (T c ) is the saturated air specific humidity, r Cm 
the unstressed canopy resistance. r cm may vary diurnally in 
a canopy depending on the atmospheric conditions but not 
the soil moisture content. However, in this study, the uns- 
tressed" transpiration is defined when the minimum uns- 
tressed canopy resistance is used in Eq. (3*3*16), which is 
a function of plant species only (from now on in the text 

r c “ r cra^ 

Alfano (1981), using data (Bruce et al., 1977, Jackson 
et al. , 1971 and 1975) for the actual evaoporation E from 

bare soil, fitted a relation between the scaling factor ^ 
and soil moisture content. The relation used in this study 
later is 

Pb =E/E p =(8 avf /Cfe) C 'f (5-5.17) 

where C g _ 50 avn to 6.50 avm and C ? ia around 5-0. 
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For a canopy, the actual transpiration E is related to 

• C 

the unstressed transpiration E u (Hanks and Ashcroft, 1980) 
through the scaling factor by 

Pc* V E u (3. 3. 18) 

should be estimated from the measured E c and E u « sever- 
al relations have been proposed for it (Hanks and Ashcroft, 
1980). Alfano (1981) proposed a relation for determining 
E c as a function of soil moisture in the root zone by using 
the data given by Denmead and Shaw (1962). Since E is sum 
of the uptakes U 1 and U^, the effective scaling factor 
is 

Pc'V'VWV- 


= Po«* Rd 1 + 

fc2*( 1-Rd 1 ) 

(3*3.19) 

where for layer i 



Pci=WV R di=<0i- 

'®pwp ^ ^ 9 t " 9 pwp ^ 

9 pwp <9 i <6 t 

or * P c i=° 

9 P wp >9 i 


or Pei*' 

9 t <9 i 

(3.3.20) 


The threshold value 0^ depends on the soil moisture and 
the value of unstressed transpiration. The relation for 0 t 
is 

®t = Vp+^fo^pwp’t^Eu^efc-Spwp’J 09 

where 0^=2520 sec/cm, C^= 0.8 and the dimension of E u is 
cm/sec . 
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/ 


3.4 COUPLING WITH PBL AND DETERMINING THE INTERFACE 

FLUTES ! 

3 . 4 .I Energy balance — determination of heat fluxes into 
soil and heat 'sTorage increment Tn canopy 


In Section 3*2, the heat fluxes into soil G„ and G v and 

c b 

the heat storage increment rate remain to be determined. 
They depend on three sub-system heat balances in the grid. 


The heat flux G^ into the bare soil surface is obtained 
from the energy balance between the bare soil surface and 
the air: 


G b = R nb.- LE b“ H b 


(3.4.1) 


where R^ is the net radiation at the surface, LE^ the la- 
tent heat flux from the surface into the air and sensi- 
ble heat from the surface to the air. The R^ is calculat- 


ed by 


< R nb- S v<’- a lb> +S L- tfT S b 4 (3*4*2) 

where S w is the solar beam from air, a^ the albedo of the 

bare soil surface, the downward long wave radiation from 
the air and 6 Stefan-Boltzmann 1 s constant. 


The heat storage increment rate in the bulk canopy layer 
Xjyj is dependent on the heat exchange between the air and 
the canopy layer as well as the canopy layer and the vege- 
tated ground surface and is expressed by 

X M= R nc -lE c ('-C l )-S w EXP(-IAI.X (; )(t-a l0 )-H 0 +H 0t , 

« (T gc 4 -T c 4) (3*4*3) 


+ 



R nc is the net radiation at the top of the canopy and is 
calculated 

R no= S w< , - a lc> tS L- 6T c 4 <5.4.4) 

where a^ c is the albedo value of the canopy. 

LE is the total latent heat from the canopy layer, 
which is contributed mainly from the canopy layer itself 
and a small portion of which come from the vegetated ground 
surface, and is the ratio of the latent heat from the 

ground surface to the total latent heat from the canopy to 
the air, which is not needed in this model with the subgrid 
parameterization. The term S W EXP(-LAI*X C ) ( 1-a lc ) repre- 
sents the part of solar radiation arriving at the ground 
surface. H is the sensible heat flux from the canopy lay- 
er and H cb is the sensible heat flux to the canopy layer 

from the vegetated ground surface. The last term on the 

/ 

right hand side of (5*4 -5) is long wave radiation exchange 
between the canopy layer and the ground surface. The heat 
flux G c into the vegetated ground surface under the canopy 
is 

G c =-LE c c 1+ S w EXP(-LAI*X 0 )(!-a lc ) -H Qb 

+ (SIT 4 -? 4 ) . (3.4-5) 

In Eqs. (5.4.1), (5-4.5) and (5.4-5), the fluxes E b , 
E c , H^, H c and H c ^ remain to be determined. The next sub- 
section will discuss the estimation for the fluxes. 
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3»4 .2 Coupling with The PBL — fluxes of momentum, heat 
and moisture from interface . 

The determination of turbulent fluxes across the atmos- 
phere-land interface is the most familiar and most widely 
dispussed aspect of the PBL parameterization problem, and 
in recent years promising new theories have improved the 
understanding (Deardorff, 1972; Randall, 1983)* In the 
GHM, Deardorff's theory is used to estimate the fluxes. It 
is a kind of similarity theory, in that the fluxes across 
the interface are related to the bulk property of the PBL. 
Now the problem is posed in the following way: if the mean 

atmospheric field such as the wind velocity u m , potential 
temperature S m and specific humidity q m within the PBL and 
the thickness of the layer h were given from GCM output, 
what should be the momentum, heat and moisture fluxes 

across <the interface when the necessary conditions on the 
land surface have been prescribed? 

Before answering this question, additional notation is 
needed to supplement those listed in 'NOTATION*. 

Y log 10 (-RI b )-3.5 

Subscripts 

a ai r 

av grid averaged variables (briefly named as 'grid') 

r c over canopy 

s the interface or b over bare soil 

lower boundary of PBL l av about grid average 

m mean value within PBL 



94 

N neutral stability 

v virtual potential temperature 

w water 

* friction velocity 
Superscripts 

* saturated condition. 

The potential temperature S is a temperature that a parcel 
of air has when it moves ad iabatically from the level with 
pressure P to the level with P=1000 mb. The virtual poten- 
t tial temperature T v is approximately defined as 

T v = S(1 .+ 0. 61 q) (3-4.6) 

The kinematic vertical heat flux F is related to the kine- 
matic sensible heat flux F,p and the kinamatic moisture flux 

p w 

f - jy0.61S F (3.4.7) 

/ 

Due to the effect of buoyancy, there are three states 
of thermal stability of the PBL: stable, unstable and neu- 

tral. The state of thermal stability of the PBL is desig- 
nated by 


stable. 

if 

(T VS - T vm> <0 

with 

F<0 

unstable , 

if 

< T vs- T vm )>0 

with 

F>0 

neutral , 

if 

< T vs- T vm ) =° 

with 

F=0 ; 


where 

T vs - S 3 (t. + 0.61q 3 ) 

(3-4. 8a) 


T v„ = 

(3.4.8b) 



or may either he 0 b or q b at the bare soil surface or 


or % * f ' the canopy layer surface, and S R is 


% = T gb {l 000/P ) n * 2RP 


or 


? c = T C (1000/P) 


0.2 PR 


(3.4.Rc) 

(3-4.8d) 


The Richardson number, which is related to the exchange 
coefficient in these three cases, is defined as 

BI b = <' h C'vm- T vs’/< T vr, u m ? ' 

The stable, unstable and neutral cases correspond to RI^ 
greater, less than or equal to zero, respectively. 


According to Deardorff's work -(1972), the friction coef- 
ficient C u and heat transfer coefficient C T for these three 
cases can be expressed as a function of RI b . For the neu- 
tral case, C uN and C„ N are given by 


/ C uN = f^“ 1 ln(O.02Sh/z o )+8.4)" 1 
C Tn = (k" 1 Rln(0.02Ph/z o )+7.3)'" 1 
For the stable case, the C u and C^, are 

C u = C uM (l “ RI b/ RI c> 

C T = C TN (l-RI b /RI c ) 


(3.4.9a) 

(3.4.9b) 

(3.4. 10a) 
(3.4.10b) 


where RI c is 3.03 and RI' b should be less than 0.9RI o 


For the unstable case, the coefficients are 


C u = ^ C uN 


-1 


2\ \-1 


- 0.23exp(0.26Y-0.03Y <; ) ) 


C T = ^ C TN 1+C u 1 “ C uN 1 ^ 


(3.4.11a) 

(3.4.11b) 


The limiting case of free convection must be considered 

For RI b in the vicinity 


under the strongly unstable case. 
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of, or exceeding, the value indicative of the commencement 
of free convection, and C q, are not known with any accu- 

racy. However, the main wind must be so small in this case 
that some inaccuracy in C u may be tolerated. The C u and 
are maintained constant for lying within the free con- 

vection regime. The regime can be identified by numerically 
testing whether C^ -1 and are less than 0.6 and 0.3 of 

their respective neutral value. For the kinematic heat 
flux, however, the above procedure should be further con- 
strained with the following condition (Townsend, 1964): 

F - > 0.19(T vs -T vm ) 4/3 (3-4.12) 


free conv. 


That means, when C u ^<0.5C u jj ^ and C,p 1 <0.?C T ^ ^ is in free 


Otherwise , 
F 


free conv. 


=0. i 9(?vs“^vti)) 4 ^ (cra*°K/s) (3-4. 13c) 


The kinematic vertical heat flux F can be partitioned 
into the kinematic sensible heat and moisture fluxes by 

(3.4.14a) 


F w -E/(f a )= (q s -q m )/(T vs -T vm )P 
'CpMS a -Sj/ (T -T )F 


(5.4.14b) 


and u # can be obtained by 


u*= u C„ 
mu 


(3.4.14c) 



G7 


Eqs. ("5. 4. 6) through (5.4.14) are applied to the bare 
soil and the bulk canopy layer, respectively. Under the 
same mean PBL conditions, different surface conditions of 
either bare soil or canopy layer will correspond to differ- 
ent sensible heat and latent heat fluxes. 


5.5 GRID AVERAGE VARIABLES 


In the previous section, the calculation of the depen- 
dent variables at the different parts in a grid was de- 
scribed. However, in order to be consistent with a GCM, 
The GHM should supply the GCM with a grid average sensible 
heat, a grid average moisture flux, a grid average surface 
temperature and a grid average friction velocity (later 
each variable preceded with grid indicates the variable are 
grid-averaged.) The basic concept developed in the GHM to 

get the grid surface temperature T is such that the calcu- 

6 

lated grid friction velocity u #av , grid kinematic vertical 
fluxes F_ . grid kinematic sensible heat and moisture flux- 
es from the ground surface, based on the grid surface temp- 
erature being calculated below and the same conditions in 
the PBL, are equal to the proportional summation of the 
corresponding values of the canopy part and bare soil part. 
Therefore, one has 


0 0 ^ A 

u *av ' 5> u *b + o u *c 

(%5. la) 

p »v * «b p b + 6 c P c 

C5.S.1b) 

P Tav = 6 b p Tb + S e P Tc 

(3.5.1c) 

P wav * V-b + ^wo 

(3.5. Id) 





QR 

p y usin/’ the formulae in Section 3-A* the grid ground sur- 
face temperature can be calculated inversely with given 

■ p 

P av’ P Tav’ P wav ’ u *av" and the Same mean PBL conditions. 


(liven F_ . the stability can be determined from 

V 


stable 

if 

P av <0 

unstable 

if 

P av >0 

neutral 

' if 

F av=° 


In each case, the grid friction coefficient may be directly 
calculated by 


C u = u -av /u n (3 - 5 - 2) 

According to Fq. (3-4-9) ♦ grid surface roughness Z Qav 
may be calculated for the neutral case as 

lnz oav = ln(0.025h)-k(1/C uN -8.4) (3-5-3) 

Grid average heat coefficient can be calculated by in- 

serting z Qav from Eq. (3-5-3) into Eq . (3.4-9t>). Then, us- 
ing Eq . (3.4.14b) and the definition of F, the grid poten- 

tial temperature at the ground surface S for the neutral 
case is 


P av = P Tav /(u 


*av C TN^ +S m 


(3-5-4) 

For the stable case, an explicit solution can be obtained 
through mathematical treatment. First, an abbreviated nota- 
tion for convenience is introduced as follows: 
y=k _1 ln(0.025h/Z oav ) 


x ~ ^ 1 vm T vav' 
u=eh/(T vmU|n z ) 


( 3 . 5 . 5 a) 

( 3 . 5 . 5 b) 

(3-5- 5c) 
(3- 5- 5d) 


a^.^-P^R 
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/ 


Prom Eqs. (^.4. 10a) and (^.4.9a), one has 
-1 


C u C uN 

n -1 


= 0 *- RI b /Rl c ) 


-UN =(y +R *4) 

C u (y+«.4)=(1 .-RI h /RI e ) 

Also from Eqs. (3*4. 10b) and (3* 4 -9b), one has 
^T^TN = 0 •” R ^ b /RI c ) 

C TN~ 1=(Ry+7 * 3) 

C T (Ry+7.3)*(l—RI b /RI c ) 

Then, eliminating the variable y from Eqs 
(3. 5* 7c) , one has 

C T (R(1.-RI b /RI c )/C u + 7.3 - 8 . 4 R) 

= <1- R V RI c> 


(3.5.6a) 
(3.5.6b) 
(3.5.6c) 

(3.5.7a) 
(3.5.7b) 
(3.5.7c) 
(3.5.6c) and 


Substituting Rl^sbx and 
( 3 . 5 .R), one has 


-1 


= u. 


av 


x/(-P av ) 


( 3 . 5 . 8 ) 

into Eq. 


(R/C ut a)-R/(C u RI c )bx = -u, av x/P ay 

+u .av bxZ/(P av BI o> <3-5-9) 

After simplifying the terms', Eq. ( 3 . 5 . 9 ) becomes: > 

u. av bx 2 /(P av RI 0 ) + lRb/(O u RI 0 )- 

U »av /P av) x -(R/C u + a) * 0 (3.5-10) 

The above equation can be rewritten as 

Ax 2 + Bx + C = 0 (3.5.11a) 

The solution of x is 

x= t-B-(B ? -4AC) 1 / 2 J/(2A) ; (3-5-llb) 

The reason for the negative sign before the square root is 
that x approaches to zero when F av goes to zero. 
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f 


The final solution should include the limitation of the 
stable case 

RI b = bx < 0 .qRI c ( 3 . 5 . 12 ) 

If the newly calculated RI b , based on the x from Eq . 
( 3 . 5 . 11 b), is greater than 0 . 9 RI C> the solution for x 
should be changed to 


x= 0.?RI c /b (3.5.13) 

Then, 

c uN = C u /( 1 •- bx/RI c ) (3.5.14) 

and z Qav can be calculated with the C u ^ from Eq. (3* 5* 3) 
Using Eq. (3.4.14b), the grid surface potential temperature 
$ av becomes 


\v ^m R Tav^ R *av x 


and T g = s av ( P/1 000. ) (3-5. 15) 

For the unstable case, the explicit solution can not be ob- 
tained v However, it is possible to seek an equation for x 
which is easy to iterate and fast to converge. Again, using 
abbreviation 




Y= Y (x) = log 10 (-RI b )- 3.5 

(3.5.16a) 



X= X (x) = exp(0.26Y-0.03Y 2 ) 

( 3 . 5 . 16 b) 

From 

Fqs 

. ( 3 * 4 . 11 a, b), one has 




C - 1 = C m _ 1 -25X 
u uN 

(3.5.17a) 



Cm" 1 = C tn " 1 -25X 

(3.5.17b) 

From 

Eqs 

. (3.4.9a,b), one has 




c uN~ 1= y +8 ‘ 4 

(3.5.18a) 



^TN = ^Y + 7 * 3 

(3.5.18b) 



101 

Substituting Eqs. (?.s.ipa,b) into Eqs. (3.5.17a,b) and 
eliminating y, one has 

■ u *av x/F av s “ 2S( 1 _R ) x+0 - *( 1 )-1 • 1 +C U " 1 R (3- 5. 1 9a) 

or aT = F(at) (3.3.19b) 

where a?= -x and F(-x) is equal to the right-hand-side of 
Fq. (3.3.19a) divided by u* av /F av . Eq . (3-5. 19b) is nonli- 
near and can only be solved by an iteration method. 

Examining Eo . (3.5.16a), it can be shown that Y will in- 
crease with increasing AT. Also, log 10 (-RI b ) is always less 
than 5*5 and thus Y is always negative. X will increase 
with increasing Y because 

dX/dY= exp(0.26Y-0.03Y 2 )(0.26-0.06Y)>0 

Since 

d F/ d (a T ) = -25(l-R)dX/dYdY/d(^T) 

and 

/ 

dY/d(AT) = -b/RI b >0 

dF/d(z*T)<0 and thus F(-^T) decreases with increasing T. 
Under the unstable state, a t is greater than zero. So, F(a 
T) is greater than zero at ^T=0. Figure 3«5.1a shows the 
scheme for the curves of the left-hand-side and right-hand- 
side in Eq. (3-9. 19b) 
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Figure 3 . 5 -la Iteration scheme 


If a simple and direct iteration is used, choose 

AT n+1 =F(4T n ) (n=1 ,2,3,4 ) and the equation may not 

converge to the true solution a Tq or may converge to it 
slowly (see Figure 3 . 5 -la). In this study, ^T n+ 1 =0-5 

(aI n +F(&T n )) (n=1,2,3 ) was chosen and the iteration did 

converge to the real solution more rapidly (see Figure 
3-5-lb) . 

After the solutionaT 0 is obtained, S av and Z Qav can be 
solved as in the stable case. Then, grid surface tempera- 
ture T can be calculated from Eq. (3*5*15) • 

O 




If the lower boundary of the PBL is located in the 
neighborhood above the roughness height of bare soil or 
canopy, the values of atmospheric conditions at this bound- 
ary may not be the same as they would be at the physical 
surface of the bare soil or canopy. In general, one may as- 
sume the air temperature at the lower boundary equal to the 
temperature of the bare soil surface or canopy. However, 
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this assumption may not be true for the specification of 
the specific humidity. This section will be devoted to for- 
mulation of the specific humidity. 

It is very difficult to directly determine the air spe- 
cific humidity q in the neighborhood above land surface. 
One of the currently-used approaches is to relate the q 

o 

with q(T ) by the scaling factor. 

3 

In terms of resistances, the potential evaporation E , 
and the actual evaporation from the' bare soil can be 
expressed by 

V * y a ( 1b*- , im )/r v (3.6.1a) 

V ? a ( ‘»tT‘lm )/r v (3.6.1b) 

where r y is the air resistance for the evaporation from the 
bare soil. The air specific humidity at the lower boundary 
of the ^BL q^ can be estimated from Eq. (3.6.1) using the 
scaling factor derived in Eq. (3*3.17) 

Pb%* + if %* >q m 

or q b =q b* if q b #<q m (3-6.2) 

For a canopy, the maximum unstressed transpiration E u , 
according to its definition, is equal to 

E u - + r c> .. (3.6.3a) 

V ? a (< «o- <, m )/r v (3.6.3b) 

where r is the air resistance above canopy and r is the 
minimum unstressed canopy resistance. Defining 
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f c = (1 + r c /r v )_1 (3.6.4) 

and using the same derivation and assumption for over 

hare soil, one has 

q o= f c l 3 c <1 c* +( ’- f c Pc^m lf 9 c *>') n 
or q 0 =q 0 * If q c ’<q m (5-6.5) 

In Eqs . ( 3 . 6.1 ) and (3*6.3), r y depends on the air specific 
humidity and temperature near the land surface. Strictly 
speaking, Eqs. ( 3 . 6 . 1 ) and (3*6.3) should be solved itera- 
tively through the PBL. However, in this study, r v is cal- 
f culated assuming the air specific humidity near the land 

surface is under saturated conditions at the land surface 


temperature . 



Chapter IV 

COMPUTATIONAL SCHEME 


4.1 METHOD FOR SOLVING SOIL MOISTURE AND TEMPERATURE RATE 

EQUATIONS 

In the GHM, a system of six basic differential equa- 
tions should be solved simultaneously grid by grid. Since 
F the six time-dependent variables can be classified into two 

basic types: fast variation (T ^ and T c ) and slow variation 
(T, -g C * and @2^* a l3ackwar ^ implicit finite difference 
method and an explicit finite difference method are applied 
to the equations with fast variation and those with slow 

change, respectively. 

/ 

The explicit approach is straight forward and discussed 

here briefly. For example, if one has nonlinear equations 

for an unknown vector A. ( j=1 , 2. . . ,m) : 

I 

dAj/dt = F ^ ( A ^ , ..., Aj^ ), 3 = 1 , . . ,m (4*1*1) 

The basic formula for solving the nonlinear system with 
the explicit finite difference scheme is 

A j n+1 = A j n+DtP j( A i n ’ A 2 n A m n ) (4.1.2) 

where n indicates the nth time step and Dt is the time in- 
terval, 1 POO seconds used in this study. 
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For the backward implicit finite difference, the dynamic 
system defined by (4.1.1) can be expressed by 

{A j n+1_A j n)/Dt = p j (An) + 2 dP j n/d V A k n+1 - A k n > 

K ( j = 1 , 2, . . .ra) (4-1.3) 

Tn the OHM, a simpler scheme is used , that is the second 
term in the right-hand-side of Eq. (4*1.3) involves only 
the fast variable. 


4.2 IMPLICIT SCHEME FOR BARE SOIL SURFACE TEMPERATURE 

The governing equation for bare soil surface temperature 
is 

dT^/dt = [2G b /(Cd 1 )-2’7r(T gb -T)J/oi = F. (3-2.1) 
Thus, the backward implicit scheme with an expansion term 
of V will produce 

(T gb n+1 -T gb n )/Dt=2G b n/ (oL c d"'i) -2 IT /( X* ) (T gb °- 



T n ) +2 / ( oL Cd 1 ) d G b n / d T gb ( T gb n+1 -T gb n ) 
-2K/(r°0(T gb n+1 -T gb n ) 

(4.2.2a) 

Then , 

T gb n+1 = T gb n +F n (....)/(l/Dt-dF/dT gb n ) 

(4.2.2b) 

where 

dF/dT gb n = 2/(dC<f, )dC b /dT gb n -2 7T/(7:cq 

(4.2.2c) 

Since 

V R nb- L V H b 



R nb= S w( , - a lb )+ V S ' T B b 4 


so , 

aC b /dT gb * dR nb /dT gb -MV' T gb - dH b /dT gb 



(4.2.3) 


and 
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dR nb /dT gb 


(4*2.4) 

Since 

H b = •?a CpC Tb u *b^b _S m^ 
LE b = -Pa LC Tb u *b^ q b“ q iJ 



and 

U *b “ ^ub U ?n 


(4.2.5) 

thus , 

dI V dT gb = ^a Cpu m ,C ub C Tb 

d V dT g b +( V s m > 



^ C ub dC Tb/ dT gb +C 

Tb dC ub^ dT gb* ' 

(4.2.6a) 

and 

Ld V dT gb * V^- C ub C Tb 

d 1b /dT gb +(< !b-'im > 

- 


*( C ub dC Tb /dT gb +C 

Tb dC W dT gb' ‘ 

(4.2.6b) 

a v /d V" d V dT g b’ dt W d V and 

d ^ub^ d ^gb are der * ved one 

by one 

as follows. 



Since S 

b ■ T gh (l000/P)°- 2Ra , 



thus , 

dS b /dT gb = (1000/P) 0,2RB 


(4.2.7) 

Since 

q b = f b q b + ^ 1 “ P b^ q m 

if 


or < 

q b=<lb* 


(4.2.8) , 

thus , 

dq b /dT gb = V q b* /dT gb = 

541 a ? b q„*/T gb 2 




lf \ ><I ” 


or 

dq b /dT gb = 0 

if 1b‘ <q m 

(4.2.9) 

where d 

^b/dTgb = 0 because P b is 

only a function 

of e t • 

dC Tb /dT 

gb and dC ub /dTg b depend on 

the stability condition. 

For the 

neutral case. 




d W dT gb= d W dT «b =0 


(4.2.10) 

For the 

stable case, C ub and C Tb 

are given by Eq. 

(3-4.10). 

So , 





d W d V = ' C '.N /RI c dR V dT gb 
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(4.2.11a) 


or =0 if Rl, =0. 9RI_ 

b c 

dC Tb/<'V = - C M /RI c' iR1 ) 1 / ' , V 

or =0 If RI b =0.°RI c (4.2.11b) 

where dRT b /dT gb * -d? vb /dT gb h g /(T vm u m 2 ) (4.2.12) 

an4 d T vb / dT gb = 0. 61 dq b /dT gb S b +d ^b/ d ^gb^ * + ®* ^ 

(4.2. 13) 

For the unstable case, there are two subcases: the normal 

case and the free convection. Under the normal case, the 
equations for dC ub /dT b an ^ dC,p b /d?g b are : 

dC ub /dT gb = C ub 2 (2SX(Y)(0.26-0.06Y)(dRI b /d'P gb )/ 

(RI b ln10) (4.2.14) 

end dC^/dT^ = (C Tb /C ub ) 2 dC ub /dT gb (4.2.15) 

Under the free convection, there are also another two sub- 
cases : 

1 ), If C^-So.SC^- 1 , then 

dC Tb /dT gb = 0 (4.2.16) 

Similarly, if C ub “^ =0. , " then 

dC ub /d,n gb = 0 • (4.2.17) 

2) if C<pb = 0* 1 ^^ v b“"^vm^ ^^ u *b’ ^ en 

dn 'tb'' aT gb= 0 - 1q(T vb- T vn, ) ' 2/, '' ( ’ u .b )aT vb/ aT gb 

(4.2.18) 
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4-3 IMPLICIT SCHEME FOR CANOPY TEMPERATURE 


The equation for T c is 

C 0 W' ,t = X M 


(3.2.3) 


The backward implicit scheme with the expansion terms of 


T c and Tg C for the canopy temperature is 


C c h c (T c n+1 -’ P c n ' /D1; * *M n+d V aT "~< T --" + ’- T -- n > 


'M ™M 
+dX M /dT (T 


gc' gc 'gc 
n+1_ T n^ 


So , 


n+1 
c ~c 


T.— =T e n + (X M n + dX M /dT (!0 n (T gc n+1 - 
T gc n))/(h c Cc/Dt “ dX M /dT c n) 


( 4 . 2 . 20 ) 


where 


dX M / dT =4 6" T %dH . /dT 
M' gc ^ gc cb' gc 


(4.2.21 ) 

The H cb can be approximately estimated by (Monteith, 1975) 

(4.2.22) 


H 


cb ?a. C v( T g C -T c^ r cb 


where is the air resistance for sensible heat between 

the canopy and the vegetated ground surface, and is approx- 
imated by 

r ch * '/(Vcb* (4.2.23) 

where C^ is the average drag coefficient inside the canopy 

and is the mean wind velocity inside the canopy. 

The drag coefficient and the average wind velocity in 

the canopy are both approximated by (Monteith, 1975) 

2 


C D =0.2(0.37/ln(h c /Z oc ))‘ 


u cb = °* 4u top 


(4.2.24) 

(4.2.25) 



is the velocity at the top of canopy and equal 


where n. 

top 

t.o 


u 


tOD 


=u* c /O.V71n(0.3h c /Z oc ) 


(4.2.26) 


In regard to dX M /dT c , it is represented by 

c =- 1 (’-C 1 )dE c /aT t ,-dH c /dT c+ dH ob /dT c -86'T o :5 

(4.2.27) 


d E / dT and dH /dT 
c c c 

d? l/ dT gb and dH^/dT 
T gb’ C Tb ' C ub ’ u *b ’ 

d W dT gb’ d W dT 

d0 -b ' /|,T ab’ Pb* ? Tb 
"o' T C r ^c ' C uc' u 

dC „c/' IT C ' d W d V 

f c Pc* F Tc and F wc- 


has exactly the same expressions as 
provided that the quantities , R^, 

S h» q b’ dq b/ dT gb’ dS b/ dT gb' dC Tb /,dT gb’ 

gb’ dC TN/ dT gb' dRI b //dT gb’ dT vb/ dT gb’ 
and F wb are successively replaced by E Q , 

* c * S c’ q c* d ^c^ d 'c f d V d V dc Tc //dT c’ 
d W d V dRI b /dT c , d? vc /dT c , dq c */dT c , 

and d(f 0 )/d? is assumed to be zero. 



Chapter V 

RESULT AND DISCUSSION 


5.1 SIMULATION AND SENSITIVITY STUDIES 

Simulation and sensitivity studies of the GHM have been 
conducted in a non-interactive mode. Atmospheric conditions 
are generated by the NASA/GLAS GCM (Halem et al. f 1979). 
This version of the GCM has a resident GHM which has been 
used by Mintz and Sarafini (1982) to study the global cli- 
matic soil moisture distribution. The output land surface 
data of a global simulation run for 45 days between July 10 
to August 25, 1.975 were used to drive the GHM. Eight re- 

gions across the North America continent were selected for 
the simulation and sensitivity studies (see Figure 5*1)* 
The emphasis is placed on the response of the GHM to the 
atmospheric forcing in a relatively short time period under 
a variety of climatic and land surface conditions. It is 
obvious that there are no feedback effects in this non-in- 
teractive mode. However, the interactively simulated hydro- 
• * « 
logic results of the resident GHM are available for this 

study. The basic properties for these regions are listed 
in Table 5-1 * 1 - 


1 12 
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Atmospheric conditions provided by the GCM were given in 
the neighborhood of the land surface. Alfano (1981) applied 
these surface data at the anemometer height to drive the 
GHM in his study. In this study, the surface data were con- 
verted to the mean conditions in the PBL using the Dear- 
dorrf (1972) formulation. The detailed derivation to ob- 
tain the mean conditions is discussed in Appendix. 

Since the parameterization of albedo in the GHM is dif- 
ferent from that of the GCM, the absorbed net radiation in 
the GHM is different from that provided by the GCM and 
should be adjusted prior to running the model. The adjust- 
ment was discussed by Alfano (1981) and will be described 
briefly in Appendix. 

For evaluating the performance of the GHM, the results 
of grtd surface temperature, grid available moisture frac- 
tion in the root zone, grid sensible and latent heat fluxes 
obtained from this model are compared with the correspond- 
ing results furnished by the GCM. 

The analysis of the results of the sensitivity study are 
presented in terms of the daily averaged soil moistures, 
temperatures, and grid latent heat and sensible heat flux- 


es . 



11 ** 
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o 
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Figure 5. 1 Eight regions across the North America continent 
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®s 

H 

0 

pwp 

Vegetation 

type 

Soil 

type 

Bare soil 
density 

32 

19 

-475 

• 320 

H 

seasonal 

forest(7) 

silty clay 
loam (E) 

25 

33 

13 

-465 

.162 

.061 

wood 

land (3) 

sandy 
loam (B) 

25 

33 

15 

.465 

.162 

.061 

desert 

(6) 

sandy 
loam (B) 

75 

33 

17 

• 503 

.270 

.098 

grass 
land (5) 

loam 

(D) 

50 

33 

19 

-394 

.220 

. 100 

seasonal 
forest (7) 

loamy 
sand (C) 

25 

33 

21 

/ 

-394 

■ 

.100 

seasonal 

forest(7) 

loamy 
sand (C) 

25 

34 

22 

-394 

.220 

.100 

seasonal 

forest(7) 

loamy 
sand (C) 

25 „ 

35 

17 

• 503 

.270 

.098- 

‘“grass 
land (5) 

loam 

(D) 

50 


Table 5-1*1 Properties of the regions investigated 
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5*2 COMPARISON BETWEEN THE RESULTS OP THE GHM AND THE GCM 

The resident hydrologic model in the GCM (Mintz and Sa- 
rafini, 1982) differs from the GHM in many respects. The 
difference between these two models will be explained to 
aid in understanding the subsequent comparisons. 

5.2.1 Difference betwe.en the GHM and the hydrologic model 
Tn tKe" 'GCM 

In Chapters 2 and 5* the formulation of the GHM have 
,/ been described in detail. The most important differences 

between the GHM and the GCM’s resident hydrologic model are 
the implementation of features in the GHM: characteriza- 

tion of soil and vegetation, subgrid parameterization of 
vegetation density and formulation of moisture and heat 

movements in the soil layers. 

/ 

In the resident hydrologic model in the GCM, all grids 
on the land surface are considered to be identical, con- 
sisting of a uniform soil layer with the same maximum 
available soil moisture, hydraulic and thermal properties. 
The effects of canopy are not considered and albedo is pre- 
scribed for each grid. For the moisture movement, the wa- 
ter budget in the GCM in entire root zone is calculated in 
terms of precipitation, evaporation, runout and change of 
available moisture content. Only the available moisture 
fraction 6 avfs 


in the root zone is predicted. The evapora- 
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tion is parameterized by using potential evaporation and a 
scaling factor which is a function of the available mois- 
ture fraction. For the heat transport, ground temperature 
in the GCM accounts for average temperature within a bulk 
layer equal to the penetration of the diurnal temperature 
wave. Its change is only forced by surface air conditions. 

5.2.2 Comparing the simulation results 

The eight regions shown in Figure 5*1 are Northern Mis- 
sissippi (J=32, 1=19), Northern California (J=33, 1=13) » 

Utah (J=33» 1=15), Kansas (J=33» 1=17), Eastern Missouri 

(J=33, 1=19), Virginias (J=33, 1=21), New York (J=34, 1=22) 
and Dakotas (J=35, 1=17). The experiment period is from 

July 11 to August 25 in 1975. The following analysis is 

based on comparing the daily average results from the GHM 

/ 

under the normal conditions with those from the GCM* : The 

parameters of the normal case are listed below: 
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J 

I 

Location 

Albedo 

6 £ 
.. ^ 

r c 

d 2 

R d1 

32 

19 

Northern 

Mississippi 

Formula 

25 

75 

90 

*25 

33 

13 

Northern 

California 

Formula 

25 

75 

90 

.25 

33 

15 

Utah 

Formula 

75 

100 

40 

.50 

33 

17 

Kansas 

Formula 

50 

100 

40 

• 50 

33 

19 

Eastern 

Missouri 

Formula 

25 

75 

90 

• 25 

33 

21 

V i rginias 

Formula 

25 

75 

90 

.25 

34 

22 

New York 

Formula 

25 

75 

90 

.25 

35 

17 

Dakotas 

Formula 

50 

100 

40 

• 50 


d^ = 10 cm for all grids 

Normal conditions for eight regions 

Figures 5*2-1 through 5*2-8 summarize the results of 
comparisons, where the curves with 'square* represent the 
results from the GHM, The curves with 'plus' represent the 
results from the GCM. In the temperature plots, the curves 
with 'asterisk' mean the air surface temperature from the 
GCM. The quantities with subscript 'd' indicate the daily 
average value. In figure 5.2-1c through 5*2-8c, it is evi- 
dent that the grid surface temperatures from these two mod- 
els follow a similar trend. However, the GHM predicts 
larger diurnal variation than the GCM does. 

In the forest regions with relatively moist soil (J=32 
and I=19» Northern Mississippi; J=33 and 1=19, Eastern Mis- 



souri; J = 33 and 1=21 , Virginias; J = 34 and 1=22, New York), 
the daily averaged grid surface temperatures T , from the 
GHM and the daily averaged ground surface temperature Tg g{ j 
from the GCM are in phase and agree to each other reason- 
ably well (see Figures 5*2-1c, 5*2-5c, 5*2-6c, 5.2-7c). The 
maximum difference T gsmax between T gd and T gsd for the four 
regions ranges from 1.7 to 3*0 °K and its average Tg g ^ over 
45 days ranges from 0.78 to 1.3 °K. The daily averaged 
surface air temperature T &d is also shown in the figures 
which is nearly equal to both of these temperatures. Prom 


Table 5*2.1 the maximum difference T between T„j and 

gamax aa 

T^ varies from 1.4 to 2.9 °K and its average T ga ^ over 45 
days from 0.56 to 0.87 °K for the GHM; the maximum differ- 
ence T gamax between T ad and T ggd varies from 1.8 to 3*1 ° K 
and its average T gab over 45 days from 0.69 to 1.02 °K for 
the GCM. 


There is one region (J=33 and 1=1 3» Northern California) 
with tall canopy but relatively dry soil. During the first 
20 days, the daily averaged grid surface temperature T^^ 
from the GHM is relatively close to the daily averaged 
ground surface temperature Tg g( j from the GCM in both magni- 
tude and trend (see Figure 5*2-2c). After.- 20 days, how- 
ever, the soil is dried out and evapotranspi ration becomes 
very small. Under this situation, the moderating function 
of latent heat becomes weakened and albedo plays a more im- 
portant role in determining the surface cover temperature. 
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Since the albedo in the GHM is larger than that in the GCM, 
the daily averaged grid surface temperature predicted by 
the GHM becomes much smaller than the daily averaged grid 
surface temperature given by the GCM. From Figure 5.2-2c, 
the temperature difference increases and eventually goes up 
to 6*3 °K. From Table 5-2.1, the averaged temperature dif- 
ference T^ gb between these two temperatures over 45 days is 
2.1 °K. These two daily averaged grid surface temperatures 
from these two models are both higher than the daily aver- 
aged air temperature. But the difference between T ^ from 
the GHM and T gd of the air is smaller than that between 
^gsd ^ rom the GCM an< * ^ad‘ From Table 5-2.1, the averaged 
temperature difference T^ a ^ between T gd and T ad over 45 
days is 1 .43 °K, which is lower than the corresponding dif- 
ference T gab between T g3(J and T &d , 3-13 °K. 


For the grids with 50# covered by short canopy (J=33 and 
1=17, Kansas; J=35 and 1=17, Dakotas), the daily averaged 
grid surface temperature from the GHM T^ d and from the GCM 
^gsd a € ree reasonably well in trend and in magnitude (see 
Table 5-2.1 and Figures 5-2-4c and 5-2-8c). From Table 
5-2.1, the absolute temperature differences T^ gb between 
these two temperatures averaged over 45 days are' 1 .72 and 
1.18 °K, respectively. The daily averaged surface air 
temperature T gd shown in the figures is also nearly equal 
to both of these temperatures. From Table 5-2.1, the abso- 
lute differences T^ ab between T^ d and T &d averaged over 45 
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days are 1.46 and 0.84 °K; the absolute differences T , 

sab 

between T and T &d averaged over 45 days are 0.88 and 
0.77 °K 

There is one short canopy region (J=33, 1=15, Utah) but 
with 75# covered by bare soil where the difference between 

T^ d and T g S( j is less than 3 °K on most days. But the larger 
difference between temperatures T ^ and Tg S(J occurs on some 
days. The maximum difference T g amax /between T gd and T g3(J is 
6.8 °K and its average T ggb over 45 days is 2.90 °K (see 
Figure 5*2-3c). The large difference of these two tempera- 
tures on some days can be explained by the following rea- 
sons. This grid is 75# covered by bare soil. The grid av- 
erage surface temperature is more dependent on the bare 
soil surface temperature variation. There are heavy rain- 
falls on these days when the larger grid surface tempera- 
/ 

ture difference occurs. The rainfall will increase the 
moisture content of the surface layer in the GHM and thus 
increase the evaporation from bare soil greatly. On the 
other hand, the rainfall does not significantly increase 
the moisture content of the entire root zone in the GCM and 
the evaporation in the GCM is affected by the heavy rain- 
fall to a smaller degree. Since the bare soil temperature 
in the GHM only accounts for a thin surface layer, the 
large increase of evaporation on these days will depress 
the surface temperature increment. Thus, larger grid temp- 
erature difference between T^ d and Tg gd 


occurs. 



I 

■ 

T 

gsmax 

T 

gam ax 

T 

samax 

19 

32 

BB 

2.9 

1 .8 

13 

33 

6.3 

2.8 

6.2 

15 

33 

6.8 

6.6 

2.8 

17 

33 

4.0 

m 

2.9 

19 

g| 

1 .7 

1 .4 

2.0 

21 

33 

am 

1.5 

2.0 

22 

34 

3.0 

2.6 

3.1 

17 

35 

3-5 

2.2 

mm 


Table 5.2.1 Maximum and averaged 
between the GCM, the GHM and air 


22 


T gsb 

T gab 

^sab 

1 .08 

0.87 

0.78 

2.09 

1 -43 

3.13 

2.90 

2.20 

1 .07 

1 .72 

1 *46 

0.88 

0.78 

0.56 

0.69 

1.19 

0.57 

0.85 

1 .30 

0.72 

1 .02 

1 .18 

0.84 

0.77 


temperature differences 

surface condition (°K) 
; • 





















































Since the parameterization of soil moisture in these two 
models are quite different, a comparison can only be made 
by introducing a grid available soil moisture fraction in 
the root zone as follows: 

®av - (6,*d,^ 2 *d 2 )/(d, + d 2 ) - 0 pwp 

®avf = ®av^®avm (5*1*1), 

which is shown in Figures 5*2-1b to 5*2-8b. The soil mois- 
tures from the GHM show a more sensitive response to rain- 
fall, which correlate with the rainfall events as shown in 
Figure 5 *2-1 a to 5*2-8a. The soil moisture from the GCM has 
an overall monotonically decreasing trend during the period 
of simulation. The difference in magnitudes of the soil 
moisture near the end of the simulation period are diffi- 
cult to assess because the characterization of soil and 
vegetation type varies from cell to cell as compared with a 
unifonrf maximum available soil moisture and field capacity 
prescribed in the GCM. The results of the two-layer param- 
eterization can not be examined by using grid available 
soil moisture fraction as defined in Eq. (5.1.1). However, 
the effects on the ground surface temperature indirectly 
through evapotranspiration and albedo have been demonstrat- 
ed in the above discussion. For example, in the dry and the 
sparsely vegetated regions such as Northern California 
(<J=33» 1=13) and Utah (J=33» 1=15), the grid surface temp- 
eratures of the GHM are lower than those of the GCM. The 
effects of the soil layers on evapotranspiration will be 
discussed a3 follows. 
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The daily averaged grid evapotranspirations from the GHM 
and the GCM follow similar trend. However, the magnitude 
from the GHM in general is less than that from the GCM in 
most regions (see Figures 5.2.-1d,5*2-2d, 5*2-4d through 

5.2- 8d). Only one region with more bare soil cover (Utah) 
is exceptional (see Figure 5.2-3<f)* The averaged ratios of 
evapo transpiration E ta from the GHM to evaporation E S(J from 
the GCM over 45 days are equal to 0.78 (Northern Mississip- 
pi), 0.25 (Northern California) 1.15 (Utah), 0.87 (Kansas), 
0.73 (Eastern Missouri), 0.71 (Virginias), 0.70 (New York) 
and 0.82 (Dakotas). These ratios indicate that, in most re- 
gions, the evapotranspiration from the GHM is reduced by 
1 0# to 3 0# of those from the GCM. The only grid with the 
ratio greater than unity is located in southern Utah and 
northern. Arizona. In this region, the land surface, is char- 
acterised by 75# bare soil or a sparsely vegetated surface. 
It appears that the rainfall events in the simulation peri- 
od kept the soil moisture in the GHM larger than that in 
the GCM as seen in Figures 5*2-3a and b, and thus enhanced 
the evapotranspiration. 

The sensible heat flux of the GHM shown in Figures 

5.2- 1 e to 5*2-8e varies with a larger amplitude than those 
of the GCM, which can be explained by the reduction of la- 
tent heat flux as discussed above. In general, the sensi- 
ble heat flux is the result of the complex energy balance 
of the land surface. 



5-3 SENSITIVITY STUDIES - — 

A series of sensitivity study have been conducted for 
three sets of initial conditions, two lower boundary condi- 
tions and a number of physical parameters including vegeta- 
tion density, albedo, canopy resistance, root density dis- 
tribution and root zone thickness. The results of daily 
averaged soil moisture, temperature and fluxes are present- 
ed in Figures 5-3-1 through 5. 3-5 for 45 days and for two 
regions. These are New York (J=34, 1=22) with 6* = 0.75, 

tall canopy and more water supplied by precipitation and 
and Dakodas (J=35» 1=17) with = 0-50, short canopy and 

i 

relatively dry climate. The maximum difference and averaged 
difference in the 45 days period of daily averaged soil 
moisture, temperature and fluxes are summarized in Table 

5«3-2 for all .eight regions selected for this study. Abbre- 

/ 

viated notation and symbols are given at the beginning of 
the table. 

5-3-1 Sensitivity to initial condition 

Two sets of initial conditions have been chosen as vari- 
ations to the normal case for the eight regions. They are 
the wet condition as specified by both the soil moisture 0^ 
and ©2 at the field capacity and the dry condition by the 
soil moistures at the permanent wilting point. The normal 
initial conditions were obtained from the GCM run output 
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from July 10 to August 25 , 1975* Comparison is made be- 

tween two extreme cases. In the figures, the square symbol 
indicates the results with the wet initial condition and the 
triangle symbol indicates the results with the dry initial 
condition. 

No matter where a grid cell is, the effect from differ- 
ent initial soil moisture content in the surface layer only 
lasts approximately one week (see Figure 5-3-1a). After a 
week, the soil moistures in the surface layer starting from 
f quite different initial values tend to approach to each 

other and are subsequently determined mainly by the atmos- 
pheric conditions. The reason is that the surface layer is 
thin and the initial difference of the available water con- 
taining in this layer between these two extreme cases is 

not enough to sustain the different evapotranspiration from 

/ 

the surface layer for a long period 

On the other hand, the difference in daily averaged soil 
moisture content in the lower layer 0^^ was reduced gradu- 
ally throughout the experiment period but remains at the 
end. of the period for ail regions (see Figure 5 - 3— 1b). For 
the regions with sparsely vegetated cover and shallow root 
system such as Utah (J=33, 1=15), Kansas (J=33, 1—17) amd 

Dakotas (J=35» 1=17), the water uptake function from the 

lower layer is weak and the initial difference in ©2 is 
more difficult to be eliminated. For the regions with 



dense vegetation and deep root system, the reduction was 
more substantial. This behavior implies that an influence 
of initial moisture conditions in 0^ is more significant 
and lasting than in 9^. In fact, grid evapotranspiration 
(Figure 5*3-1f)» grid averaged sensible heat flux (Figure 
5*3-1g) and layer average temperature (Figure 5-3-le) also 
exhibit a large difference for the whole experiment period. 

The effect of the initial soil moisture conditions on 
bare soil surface temperature only occurs during the first 
week (see Figure 5-3-1 c) , which is consistent with the 
trend of soil moisture in the surface layer (see Figure 
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Table 

5.3.2 

Summary of sensitivity experiment results 

d 2 


thickness of the second layer 

«0 


vegetation fraction 

r c 


canopy resistance (sec/m) 

R d1 


root density distribution 

Lower 

B.C 

lower boundary condition in layer 3 

LEM 


max (60L| E td ( 1 )-E td (2)j ) (ly/min) 

LEB 


aver (60L|E td ( 1 )-E td (2)| ) 
over 45 days (ly/min) 

HSM 


max (60|H td (l )_H td (2)l ) (ly/min) 

HSB 


aver ( 60| H td ( 1 )-H td ( 2)1 ) 
over 45 days (ly/min) 

TBARM 


max ( |T d (l )-T d (2)l ) °K 

TBARB 


aver ( lT d ( 1 )-T d ( 2)1 ) 
over 45 days °K . 

TCM' 


max (lT cd (1)-T cd (2)l ) °K 

TCB 


aver ( lT cd ( 1 )-T cd ( 2)1 ) 
over 45 days °K 

TGBM 


max (lT gbd (1)-T gb<j (2)t) ° K 

TGBB 


aver ( I T gbd ( 1 )-T gbd ( 2)1 ) 
over 45 days °K 

VMC1M 


max ( < 0 ld (l )- 0 1 d ( 2 )f ) . 

VMC1B 


aver ( d (1 )-© 1 d ( 2 )^ ) over 45 days- 

VMC2M 


max ( l© 2 d (l )- 02 d^ 2 ^ ) 

VMC2B 


aver ( le 2 d (l )- 0 2 d ( 2 )| ) over 45 days 


(x) 

indicates the extreme conditions. 


aver means average. 
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(a) Values of VMC1M 


I 

J 

6 c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.0827 

0.0190 

0.0218 

0.0135 

0.0323 

0.0045 

13 

33 

0.0236 

0.0059 

0.0061 

0.0003 

0.0094 

0.0006 

15 

33 

0.0777 

0.0213 

0.0061 

0.0009 

0.0041 

0.0034 

17 

33 

0.0726 

0.0223 

0.0213 

0.0018 

0.0134 

0.0064 

19 

33 

0.0902 

0.0171 

0.0156 

0.0096 

0.0313 

0.0022 

21 

33 

0.0838 

0.01 57 

0.0175 

0.0076 

0.031 1 

0.0023 

22 

34 

0.0883 

0.01 34 

0.0126 

0.0061 

0.0170 

0.0016 

17 

35 

0.0804 

0.0186 

0.0185 

0.0027 

0.0121 

0.0028 




(b) Values of 

VMC1B 



I 

J 

6 c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.0419 

0.0062 

0.0089 

0.0039 

0.0127 

0.0020 

13 

33'' 

0.0051 

0.0010 

0.0017 

0.0001 

0.0034 

0.0001 

15 

33 

0.0273 

0.0076 

0.0021 

0.0002 

0.0017 

0.0007 

17 

33 

0.0293 

0.0079 

0.0098 

0.0007 

0.0062 

0.0013 

19 

33 

0.0503 

0.0059 

0.0072 

0.0020 

0.0112 

0.0006 

21 

33 

0.0451 

0.0062. 

0.0073 

0.0021 

0.0128 

0.0008 

22 

34 

0.0505 

0.0050 

0.0057 

0.0019 

0.0091 . 

0.0005 

17 

35 

0.0451 

0.0104 

0.0110 

0.0010 

0.0061 

0.0010 



(c) values of VMC2M 


I 

J 


Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.1005 

0.0094 

0.0494 

0.0265 

0.0262 

0.0265 

13 

33 

0.0290 

0.0004 

0.0063 

• 0.0091 

0.0045 

0.0003 

15 

33 

0.0346 

0.0021 

0.0065 

0.0075 

0.0093 

0.0032 

17 

33 

0.1220 

0.01 37 

• 0.0520 

0.0245 

0.0313 

0.0663 

19 

33 

0.1071 

0.0065 

0.0517 

0.0477 

0.0286 

0.0050 

21 

33 

0.0995 

0.0077 

0.0439 

0.0423 

0.0259 

0.0062 

22 

34 

0.1002 

0.0081 

0.0449 

0.0435 

0.0260 

0.0050 

17 

35 

0.0805 

0.0130 

0.0425 

0.0210 

0.0210 

0.0695 




(d) Values of 

VMC2B 



I 

J 

V 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

3? 

0.0621 

0.0056 

0.0286 

0.0157 

0.0152 

0.0130 

13 

33 

0.0248 

0.0002 

0.0025 

0.0027 

0.0025 

0.0001 

15 

33 

0.0281 

0.001 5 

0.0047 

0.0053 

0.0066 

0.0023 

17 

33 

0.0607 

0.0068 

0.0274 

0.0133 

0.0152 

0.0395 

19 

33 

0.0696 

0.0043 

0.0313 

0.0276 

0.0179 

0.0024 

21 

33 

0.0566 

0.0043 

0.0237 

0.031 1 

0.0140 

0.0025 

22 

34 

0.0553 

0.0044 

0.0237 

0.0282 

0.0140 

0.0020 

17 

35 

0.0464 

0.0070 

0.0261 

0.0089 

0.-0130 

0.0350 
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(e) ' 

Values of 

TGEM ( 

°K) 


I 

J 

<0 (*) 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

5.7000 

1 .9001 

1 .1001 

0.5999 

1 .5999 

0 . 2002 

13 

33 

4.7000 

2.9001 

0.8000 

0.1001 

1 .6001 

0.6001 

15 

33 

5.3999 

4 . 5000 

0.7000 

0.2000 

0 . 5000 

0.8999 

17 

33 

3-5000 

2.2000 

• 1 .0000 

0.1001 

0.7002 

0.5000 

19 

33 

4.5000 

1 .8999 

1 .4001 

0.3000 

2 . 0000 

0 . 2000 

21 

33 

4.3000 

1 .7000 

1 .1001 

0 . 5000 

1 .6001 

0.2002 

22 

34 

4 . 5000 

2.2000 

1 .3000 

0.6001 

1 .6001 

0.2000 

17 

35 

2.5000 

2.2000 

0.7000 

0.1001 

0.5000 

0.3000 




(f) Values of 

TGBB ( 

°K) 


I 

J 


Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 < 

1 .6689 

0.6067 

0 .2600 

0.0889 

0.3733 

0.0578 

13 

33 

2.5778 

1 .6400 

0.1356 

0.01 1 1 

0.2134 

0.2378 

15 

33 

1 .871 1 

1 .5667 

0.1556 

0.0244 

0.1467 

0.21 1 1 

17 

33 

1 .4778 

M 

0.9000 

0.2889 

0.0245 

0.2178 

0.1378 

19 

33 

1 .7088 

0.6200 

0.2578 

0.051 1 

0.4733 

0.0356 

21 

33 

1 .2800 

0.5200 

0.2244 

0.0422 

0.3444 

0.0400 

22 

34 

1 .6822 

0.6378 

0.2022 

0.0733 

0.2756 

0.0489 

17 

35 

1 .0778 

0.8956 

0.2533 

0.0222 

0.1800 

0.0956 


(*) 

under 6 C 

. term means the 
* 

given values indicate 


difference of grid surface temperature between afforesta- 
tion and desertification cases. 
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(g) Values of TOM (°K) 


I 

J 

6 

c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 


0.8000 

1 -5999 

0.1001 

0.3000 

0.3999 

13 

33 


0.6001 

0.5000 

0.3000 

0.6001 

0.1001 

15 

33 


2.3000 

2.3000 

1 .2000 

1 .5999 

0.6001 

17 

33 


1 .0000 

. 2.7000 

0.0000 

0.8000 

0.8999 

19 

33 


0.5000 

0.8000 

0.6001 

0.5000 

0.2002 

21 

33 


0.4001 

1.2002 

0.6001 

0 . 2000 

0.1001 

22 

34 


0.5000 

0.7000 

0.7000 

0.4001 

0.1001 

17 

35 


’ 0 . 7000 

2 . 1 001 

0.0000 

0.3000 

0.3000 




(h) Values of 

TCB (° 

K) 


I 

J 

c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 


0.2778 

0.6356 

0.0044 

0.0289 

0.0689 

13 

3^ 


0.4045 

0.0844 

0.0889 

0.0689 

0.0422 

15 

33 


0.9778 

0.8578 

0.25H 

0.3222 

0.2000 

17 

33 


0.5156 

1 .51 1 1 

0.0000 

0.0689 

0.1578 

19 

33 


0.2956 

0.4800 

0.2244 

0.0800 

0.0222 

21 

33 


0.2356 

0.4533 

0.1155 

0.0178 

0.0089 

22 

34 


0.2533 

0.4133 

0.1111 

0.0356 

0.0111 

17 

35 


0.4733 

1 .1844 

0.0000 

0.0200 

0.0333 
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(i) 

Values of 

TBARM 

(°K) 


I 

J 

Sc 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

1 .9001 

0.5000 

0.6001 

0.1001 

0.2002 

0.2002 

13 

33 

3-3000 

1 .1001 

0.1001 

0.2002 

0.1001 

0.5000 

15 

33 

5.1001 

2 . 0000 

0.1001 

0.1001 

0.1001 

0.5000 

17 

33 

5.3000 

1 .4001 

. 0.8000 

0.1001 

0.3000 

0.3999 

19 

33 

0.8000 

0.5000 

0.5000 

0.4001 

0.3000 

0.1001 

21 

33 

0.8999 

0.4001 

0.5000 

0.3000 

0.2002 

0.1001 

22 

34 

0.8000 

0.6001 

0.4001 

0.3000 

0.2002 

0.1001 

17 

35 

6.3000 

1 .4001 

0.6001 

0.1001 

0.2002 

0.3000 




(j) 

Values of 

TBARB 

(°K) 


I 

J 

«c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.9089 

0.3622 

0.3844 

0.0156 

0.1400 

0.0222 

13 

33 < 

2.351 1 

0.7089 

0.0356 

0.0844 

0.031 1 

0.2467 

15 

33 

3.7533 

1 .2556 

0.0622 

0.0222 

0.0489 

0.2000 

17 

33 

4.0222 

0.8644 

0.4733 

0.0067 

0.1222 

0.1022 

19 

33 

0.3956 

0.3667 

0.3067 

0.1444 

0.1 511 

0.01 56^ 

21 

33 

0.2533 

0.2933 

0.2933 

0.0489 

0.1 178 

0.01 1 1 

22 

34 

0.2955 

0.3311 

0.2956 

0.0667 

0.0822. 

0.0244 

17 

35 

4.4977 

0.8978 

0.3933 

0.0067 

0.1089 

0.0822 
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(k) Values of LEM (ly/min) 


I 

J 

6 

c 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.2481 

0.0336 

0.1146 

0.0145 

0.0492 

0.0394 

13 

33 

0.1516 

0.0143 

0.0612 

0.0491 

0.0220 

0.0026 

15 

33 

0.2088 

0.0676 

0.0237 

0.0084 

0.0125 

0.01 31 

17 

33 

0.2716 

0.0515 

0.0991 

0.0043 

0.0349 

0.0325 

19 

33 

0.2450 

0.0342 

. 0.1072 

0.0864 

0.0561 

0.0203 

21 

33 

0.1766 

0.0372 

0.1159 

0.1225 

0.0378 

0.0062 

22 

34 

0.2447 

0.0347 

0.1102 

0.0725 

0.0530 

0.0032 

17 

35 

0.2320 

0.0536 

0.0873 

0.0037 

0.0207 

0.0076 




(1) Values of 

LEB (ly/min) 


I 

J 


Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.0721 

0.0171 

0.0622 

0.0023 

0.0109 

0.0079 

13 

33 

0.0314 

0.0027 

0.0109 

0.0116 

0.0066 

0.0007 

15 

33’ 

0.0648 

0.0207 

0.0074 

0.0031 

0.0045 

0.0025 

17 

33 

0.0813 

0.0181 

0.0390 

0.0009 

0.0094 

0.0060 

19 

33 

0.0817 

0.0135 

0.0552 

0.0278 

0.0164 

0.0034 

21 

33 

0.0510 

0.0144 

0.0471 

0.0152 

0.0087 

0.0010 

22 

34 

0.0675 

0.0151. 

0.0479 

0.0167 

0.0091 

0.0006 

17 

35 

0.0732 

0.0171 

0.0351 

0.0008 

0.0071 . 

0.0016 
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(m) Values of HSM (ly/min) 


I 

J 

ff e 

Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.1492 

0.0993 

0.1051 

0.0070 

0.0402 

.0.0352 

13 

33 

0.1476 

0.1 106 

0.0500 

0.0431 

0.0191 

0.0148 

15 

33 

0.1400 

0.1049 

0.0176 

0.0059 

0.0079 

0.0270 

17 

33 

0.1726 

0.0809 

0.0679 

0.0026 

0.0255 

0.0308 

19 

33 

0.1965 

0.1110 

0.0965 

0.0759 

0.0526 

0.0168 

21 

33 

0.1122 

0.0882 

0.101 1 

0.1051 

0.0332 

0.0124 

22 

34 

0.0955 

0.1097 

0.1016 

0.0656 

0.0447 

0.0119 

17 

35 

0.1104 

0.0830 

0.0654 

0.0018 

0.0160 

0.01 54 




(n) Values of 

HSB (ly/min) 


I 

J 


Albedo 

r c 

d 2 

R d1 

Lower 

B.C. 

19 

32 

0.0459 

0.0411 

0.0570 

0.0012 

0.0075 

0.0069 

13 

33' 

0.0914 

0.0892 

0.0096 

0.0102 

0.0056 

0.0077 

15 

33 

0.0530 

0.0335 

0.0054 

0.0016 

0.0025 

0.0057 

17 

33 

0.0474 

0.0309 

0.0317 

0.0005 

0.0059 

0.0059 

19 

33 

0.0572 

0.0509 

0.0518 

0.0246 

0.0140 

0.0034 

21 

33 

0.0323 

0.0355 

0.0439 

0.0131 

0.0058 

0.0010 

22 

34 

0.0420 

0.0454 

0.0446 

0.0146 

0.0067 

0.0013 

17 

35 

0.0398 

0.0355 

0.0305 

0.0004 

0.0045 

0.0023 
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5*3*2 Sensitivity to vegetation density 

The large scale subgrid parameterization of vegetation 
density is expressed in terms of vegetation fraction or 
bare soil fraction which varies seasonally and spatially 
with latitude. 

Two extreme cases of land surface cover have been inves- 
tigated in addition to the normal case. One represents des- 
ertification with vegetation fraction ^ c =0 and the other 
afforestation with 6 c =1 . In Figure 5*3-2, the square sym- 
bol represents the result of afforestation and the triangle 
symbol represents the result of desertification. The normal 
case is represented by a solid line. It Is evident that 
the density of land surface cover have significant effects 
on all soil moistures, temperatures and fluxes. 

/ 

The vegetation density exerts the most important influ- 
ence on evapotranspiration and soil moisture as shown in 
Figures 5*3-2a, b and f and Tables 5.3-2(a), (b), (c), (d), 
(k) and (1). The overall evapotranspiration for the case 
of afforestation in the New York region (J=34, 1=22) is 
almost the same as the normal case, while for the deserti- 
fication case it is reduced by an average of 0.06.75 ly/min 
or approximately 40 to 50$ (Figure 5*3-2f). However, the 
variation in the magnitude of evaporation, as modulated by 
precipitation events, is substantially reduced by affores- 
tation. The maximum difference between these two cases can 
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be as large as 0.2447 ly/rain (Table 5*3-2(k)) because of 
the difference in water uptake patterns. Evaporation from 
the bare soil is practically cut off when the surface layer 
becomes dry, while the vegetated cover can continue to 
transpire by uptaking water from the root zone. The differ- 
ence in water uptake patterns leads to very different soil 
moisture profiles in the two soil layers for these two ex- 
treme cases as shown in Figures 5-3-2a and b. The results 
are similar for the rest of the regions as indicated in Ta- 
bles 5-3-2(k), (1) except the Northern California region 
where the level of soil moisture is low during this period 
(see Figure 5.2-2b). As a result, the maximum difference of 
the evapotranspi ration in Northern California is 30# less 
and the average difference is 55# less than those of the 
New York region. The difference in soil moisture exhibits 
the same trend as listed in~Tables 5*3-2(a) and(b). 

Canopy cover has quite different thermal properties from 
bare soil. Moreover, plants take up additional water from 
the lower layer, which consumes latent heat. Therefore, the 
energy balance in these two extreme cases must be very dif- 
ferent. It results in a great changp in the sensible heat 
flux, heat flux into the soil, layer-average temperature 
and grid surface temperature. Figure 5*3-2g shows the dif- 
ference in sensible heat. In Tables 5«3-2(m) and (n), the 
maximum differences of daily averaged sensible heat in 45 
days HSM vary from 0.086 to 0.20 ly/min and the averaged 



differences of daily averaged sensible heat over 45 days 
HSB vary from 0.052 to 0.091 ly/min. Figure 5*3-2e shows 
the difference in layer-average temperature. From Tables 
5 • 3—2 ( i ) and (j), the maximum difference of daily averaged 
layer- average temperature in 45 days TBARM varies from 0.8 
to 5*3 °K and its average over 45 days TBARB varies from is 
0-25 to 4-5 °K for all regions. Figure 5«3-2c shows the 
difference in the grid surface temperatures for these two 
cases. In Tables 5*3-2(e) and (f)TGBM ranges from 2.5 to 
5.7 °K and TGBB from 1 . 1 to 2.6 °K (note TGBM or TGBB here 
describes maximum difference or averaged difference over 45 
days of daily averaged grid surface temperatures between 
two cases instead of bare soil surface temperature). The 
above values indicate that the afforestation reduces aver- 
age surface temperature about 2°K and the maximum differ- 
ence in daily averaged surface temperature can be two to 
three times the above values. It should be noted that the 
subgrid change of vegetation density affects all the param- 
eters and hence outputs because the canopy implies differ- 
ent root system, canopy resistance, albedo and so on when 
compared with the bare soil. In this sense, an accurate 

estimation of o is very important in 'a GHM. 

c 
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5 *3* 3 Senaitlvity to canopy resistance 

In this study, two selections of minimum canopy resis- 
tance are added to the normal case. The values of minimum 
canopy resistance can depend on the species and physiologi- 
cal conditions of plants and can vary greatly. Table 3. 3-3 
only lists the values that are used for this study. 

Canopy resistance controls transpiration through leaves 
and thus governs water uptake from root zone. Prom the re- 
sults of sensitivity study, it is found that the change in 
the canopy resistance primarily affects the transpiration 
and soil moisture in the lower layer. The following dis- 
cussion is based on comparison of the results with the two 
extreme values of canopy resistance. 

Canopy resistance has only a small effect on soil mois- 

/ 

ture content in the surface layer as shown Figure 5*3-3a. 
In Tables 5*3-2(a) and (b), VMC1M varies from 0.006 to 0.02 
and VMC1B from 0.C017 to 0.011. These are small differenc- 
es . 

The moisture content in the lower layers is influenced 
more which tends to increase with time during the experi- 
ment period. Figure 5*3-3B gives a typical trend for the 
eight regions with tall canopy or short canopy. In Tables 
5 • 3-2 (c) and (d) VMC2M varies from 0.006 to 0.052 and VMC2B 
from 0.00257 to 0.031 • 
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The effect on the temperatures is generally weak. The 
temperature of the bare soil portion in a grid is affected 
indirectly and shows little difference between the two dif- 
ferent canopy resistances (see Figure 5-3-3c). On the other 
hand, the canopy temperature is directly affected by dif- 
ferent latent heat flux. However, the difference in the 
temperature is not significant because the heat capacity of 
canopy is greater than that of the thin soil surface layer, 
especially for tall canopies. Figure 5*3-3d shows the typi- 
cal change of canopy temperature in all regions with tall 
canopy and short canopy. From Tables 5*3-2(g) and (h), TCM 
varies from 0.50 to 2.70 °K and TCB from 0.084 to 1.51 °K. 
The layer-average temperature has less than 1 °K difference 
during the whole operation period, which only slightly in- 
creases. It is not clear whether larger difference will oc- 
cur after a longer period experiment. In Tables 5-3-2(i) 
and (j), TBARM varies from 0.1 to 0.8 °K and TBARB from 
0.036 to 0.47 °K. 

The evapotranspiration is sensitive to the variation of 
the canopy resistance. Figure 5«3-3f demonstrates the sig- 
nificant difference in evapotranspiration. Sensible heat 
always acts to compensate the latent heat in the energy 
balance and a bigger difference also occurs (see Figure 
5.3-3g)* These results are typical for the regions for 
both tall canopy and short canopy. In Table 5.3-2(f), EVM 
varies from 0.024 to 0.12 ly/min and EVB from 0.007 to 


0.062 ly/min. Prom Tables 5*3-2 (m) and (n) HSM varies from 
0.018 to 0.10 ly/min and HSB from 0.005 to 0.057 ly/min. 


symbols 
in Fig. 

solid 

curves 

(normal) 

square 

triangle. 

I 

J 

r 

(s/m) 

r c 

(s/m) 

r c 

(s/m) 

19 

32 

75 

100 

200 

13 

33 

75 

100 

200 

15 

33 

100 

50 

200 

17 

33 

100 

50 

200 

19 

33 

75 

100 

200 

21 

33 

75 

100 

200 

22 

34 

75 

1f)0 

200 

17 

35 

100 

50 

200 


Table 5*3*3 Different selection of r 

c 
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5«3*4 Sensitivity to albedo 

In the GHM, the albedo varies according to the formula- 
tion given in Section 2.6. The albedo affects the energy 
balance on the land surface. In this study, albedo was 
varied by increasing the normal value by 15# and decreasing 
it by 5#* The results shown in Figure 5*3-4 and Table 
5.3-2 indicate that the sensible heat fluxes change in all 
regions. The difference between the two extreme cases can 
reach the order of o.l ly/min. HSM varies from 0.083 to 
0.11 ly/min and HSB from 0.031 to 0.089 ly/min. 

The associated changes in temperatures are as follows. 
First, the bare soil surface temperature has been affected 
to several degrees. In some regions with low level of soil 
moisture, where the evaporation can not compensate the var- 
iation / of sensible heat, the bare soil surface temperature 

differences can become large, e.g. the Utah region (1=1 5 * 

9 

J=33) as listed in Table 5-3-4(e). In Tables 5*3-2(e) and 
(f), TGBM varies from 1.7 to 4-5 °K and TGBB from 0.60 to 
1.60 °K. Second, the effect on the layer-average soil 
temperature increases with time as shown in Figure 5*3-4e. 
In Tables 5*3-2(i) and (j), TBARM varies from 0.4 to 2.0 °K 
and TBARB from 0.29 to 1.26 °K. For canopy temperature, 
the difference between the two cases is smaller than for 
bare soil surface temperature because of the larger heat 
capacity of the canopy (see Figure 5*3-4d). In Tables 
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5*3-2 (g) and (h), TCM varies from 0.4 to 2.5 °K and TCB 
from 0.24 to 0*98 °K. 

Since albedo only influences the absorbed solar radia- 
tion, the latent heat is only indirectly affected through 
the energy balance. The difference in latent heat is about 
one third of the difference in sensible heat. Figure 5-3-4f 
shows the typical results of the eight regions for grid av- 
eraged latent heat flux. Since evapotranspiration is af- 
fected weakly, the moisture contents in both layers would 
be influenced weakly too. Figure 5*3-4a and b show the 
typical change in 0^ d and e 2d . In Table 5. 3-2 (a), VMC1M 
varies from 0.006 to 0.022 and VMC1B from 0.005 to 0.01 04- 
In Table 5*3-2(b), VMC2M varies from 0.0004 to 0.014 and 
VMC2B from 0.0002 to 0.007. 

/ 

5.3.5 Sensitivity to root density distribution , root zone 
thickness and lower boundary condition 

The remain sensitivity studies were conducted for the 
distribution of root density, thickness of root zone and 
boundary condition at the bottom of root zone. Since the 
influence of these parameters are in general small, the re- 
sults are combined for comparison with the normal case (see 
Figures 5«3-5a to g) . However, the maximum differences and 
the average differences in 45 days for all regions are sum- 
marized in Tables 5*3-2(a) to (n). 
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Root density is a parameter describing the relative root 
distribution in the two layers. It partitions transpiration 
into two uptakes, one from the surface layer of 10 cm and 
the other from the lower layer. For tall forest, most of 
the roots are located in the lower layer and the root den- 
sity fraction in the surface layer R d1 for the normal case 
is designed to be equal to 0.25. For short plants or 
crops, a large part of the roots gathers in a shallow sur- 
face layer and R^ for the normal case is designed to be 
equal to 0.5* However, this distribution also changes with 
region and season. In this study, R d1 =0.5 for tall canopy 
regions and R dl =0.75 for short canopy regions are used as a 
variation of the normal case. 

As a result of the root density distribution, the water 

uptake for transpiration is partitioned, which has a direct 

/ 

influence on the amount of available soil moisture stored 
in the two layers or the root zone. However, the mean soil 
moisture contents in these two layers also depend on the 
thicknesses of these layers. In this study, the surface 
layer is fixed at 10 cm in thickness. In the normal case, 
the thickness of the lower layer is taken as. 40 cm for 
short. canopy and 90 cm for tall canopy. For .the sensitivity 
study, 90 cm for short canopy and 40 cm for tall canopy 


were used. 
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At the bottom of root zone, there are normally two kinds 
of boundary condition that are used in modeling the soil 
moisture transport: constant soil moisture content and 
zero soil moisture content gradient. In this study the 
condition of constant soil moisture content at the field 
capacity is used for the normal case and the condition of 
zero gradient is used as a variation. 

As seen in Figure 5*3-5» significant changes occur only 
in the soil moisture content in the lower layer, latent and 
sensible heat fluxes for both tall and short canopy. In the 
figures, the canopy type in a region can be found in Table 
5.1.1. 


In Figure 5«3-5b» the soil moisture in the lower layer 

becomes smaller as the root density fraction is decreased 

/ 

because of the reduction in water uptake. However, the ten- 
dency for the thickness of the lower layer is difficult to 
assess due to the fact that the soil moisture depends not 
only on the uptake but also on the change of the interfa- 
cial fluxes. The soil moisture also becomes small in the 
second layer for the zero gradient condition because an up- 
ward interfaci'al flux is not possible in this case. 

In Figure 5*3-5f and g, the effects on the latent and 
sensible fluxes are small and more for the tall canopy than 
for the short canopy. 
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5.3*6 Summary of the sensitivity study 

The above analysis has shown that different parameters 
have different effects on the water transport and heat 
transfer. Por example, canopy resistance produces a strong 
influence on evapotranspiration and soil moisture content 
but has a weaker effect on the temperatures. Conversely, 
albedo shows a greater effect on the energy balance and 
temperature but has a small influence on evapotranspiration 
and soil moisture content. The vegetation density has 
strong effects both on the water transport and heat trans- 
fer and changes temperatures and moisture contents. But, 
the parameters such as d£, and the lower boundary con- 

dition only exert minor influences on the land surface pro- 
cesses. The sensitivity of the parameters selected for this 
study a,re summarized in Table 5*3*4 in the order of their 
significant changes as listed in Table 5*3-2. The impor- 
tance of a parameter is indicated by a number. For example, 
^ has the most important effect on all variables listed in 
the first column while albedo has the least effect on 0 O . 
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Figure 5.2-id Comparison of latent heat Figure 5.2-ie Comparison of sensible heat 
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Figure 5 • 2-7d Comparison of latent heat Figure 5.2-7e Comparison of sensible heat 
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Figure 5‘3-lb Sensitivity to initial condition 




Figure 5.3-ld Sensitivity to initial condition 
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Figure 5.3-1e Sensitivity to initial condition 
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Figure 5.3-2a Sensitivity to vegetation fraction 



VEGETATION FRACTION ol VEGETATION FRACTION 

fc DESERTIFICATION a; A DESERTIFICATION 

D AFFORESTATION £ q AFFORESTAT TON 

- NORMAL - NORMAL 




Figure 5-3-2b Sensitivity to vegetation fraction 
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Figure 5.3-2e Sensitivity to vegetation fraction 
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Figure 5«3-2f Sensitivity to vegetation fraction 
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Figure 5. 3-3a Sensitivity to canopy resistance 
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Figure 5 . 3-3b Sensitivity to canopy resistance 




Figure 5.3~3 c Sensitivity to canopy resistance 
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Figure 5*3-3e Sensitivity to canopy resistance 
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Figure 5.3-3f Sensitivity to canopy resistance 
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Figure 5.3-3g Sensitivity to canopy resistance 
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Figure 5«3-^a Sensitivity to albedo 
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Figure 5.3-^b Sensitivity to albedo 
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Figure S.B-^c Sensitivity to albedo 
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Figure 5.3-^d Sensitivity to albedo 
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Figure S.S-^e Sensitivity to albedo 
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Figure 5.3-^f Sensitivity to albedo 
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Figure 5.3-^g Sensitivity to albedo 
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Figure 5.3-5a Sensitivity to lower boundary condition 
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Figure 5.3-5b Sensitivity to lower boundary condition, 
root density distribution and root zone thickness 




Figure 5.3-5c Sensitivity to lower boundary condition, 
root density distribution and root 2one thickness 
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Figure 5-3-5d Sensitivity to lower boundary condition, 
root density distribution and root zone thickness 
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Figure 5.3-5e Sensitivity to lower boundary condition, 
root density distribution and root zone thickness 
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Figure 5 v 3-5g Sensitivity to lower boundary condition, 
root density distribution and root zone thickness 






Chapter VI 

CONCLUSION AND FUTURE STUDY 


6.1 CONCLUSION 

This study is a continuation in the development of a 
large-scale ground hydrologic model at the University of 
Connecticut for U3e in the NASA-GLAS atmospheric general 
circulation model (Lin et al., 1978). The major refine- 
ments of the GHM, which include addition of a canopy layer 
in the subgrid parameterization of vegetation density and 
coupling of the GHM with the Deardorff (1972) version of 
planetary, boundary layer parameterization, have realisti- 
cally enhanced the hydrologic feature in the GHM and im- 
proved the estimate of momentum, heat and moisture exchange 
at the atmosphere- land interface. From the 45-day simula- 
tion study and the comparison with the resident GHM in the 
GCM, the following conclusions are made: 

1 ) The actual evapotranspiration predicted by the GHM is 
in general reduced by 10# to 30# from this given by the 
GCM. The resident GHM in the GCM is known to overestimate 
the actual evapotranspiration by the 'bucket' model. 
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2) The soil moistures in the two layers of root zone in 

% 

the GHM are more sensitive to the precipitation events than 
that from the resident GHM. The variation of soil moisture 
appears to correlate with the precipitation events. The 
soil moisture in the resident GHM is constantly reduced 
throughout the experiment period, which might be caused by 
the overestimate of evaporation. 

3) The results of soil moisture, temperature and fluxes 
from the GHM represent cell to cell variation as a result 
of the more realistic characterization of land surface cov- 
er in the model. 

4) The GHM is a two-layer model for soil moisture and 
heat transports. Soil moisture and temperature in the sur- 
face layer are more sensitive to the diurnal forcing of the 
atmosphere, while the lower layer provides an adequate 
storage of soil moisture and heat for seasonal changes. 

5) By coupling with the PBL in the GCM, the GHM provides 
a more reasonable exchange mechanism in the atmospheric 
boundary layer. For such a large horizontal scale, it is 
not conceivable that the atmospheric cpnditions are uniform 
over the whole grid at the 'anemometer height'. 

The sensitivity study of the selected parameters pro- 
vides the following conclusions: 
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1 ) The effect of. initial moisture condition in the sur- 
face layer lasts in the order of five days even under the 
extreme conditions of the field capacity and the permanent 
wilting point. However, the difference in the soil mois- 
ture content in the lower layer which comprises most of the 
root zone persists for all 45 days in the experiment peri- 
od. Therefore, it is extremely important for a short-term 
numerical experiment that the initial volume of root zone 
soil moisture be accurately estimated. Even for a longer 
experiment , such as seasonal or annual, it is likely that 
the initial noises in the root zone soil moisture are dif- 
ficult to eliminate because such a long term (over 45 days) 
is required. 

2) The relative importance of the parameters investigat- 
ed in this study is summarized in Table 5 . 5 - 4 . It indicates 
that the subgrid parameterization of vegetation density is 
the most significant parameter and the root density distri- 
bution, root zone thickness and lower boundary condition 
have minor effects on the response of the GHM. 

5) The vegetation density exerts influence on all soil 
moistures, temperatures and fluxes. However, the most im- 
portant influence is on evapotranspiration and soil mois- 
ture. Afforestaion tends to increase evapotranspiration to 
a more extent for a sparsely vegetated wet region than a 
dry region. 



6.2 FUTURE STUDY 


1 ) Field data is required for evaluating the model per- 
formance as well as the isolated processes parameterized in 
the model. In addition, the sensitivity study has shown the 
importance of accurate initial soil moisture condition in 
the root zone. Therefore, it is desirable to collect large- 
scale data sets for soil- moisture and ground temperature 
globally and regionally. 

2) This study was conducted in a non-interactive mode. 
The GHM should replace the resident GHM in a GCM to conduct 
the simulation and sensitivity studies interactively. 

3) This study was limited to a time period of 45 days 
due to the limited availability of GCM-generated atmospher- 
ic data. It is desirable to extend it to seasonal and in- 

/ 

terannual experiments. 

4) Only eight isolated grids in the North America conti- 
nent were selected for this study. Future numerical experi- 
ments should be extended globally or to include contiguous 
grids in large regions with different climatic conditions. 



Appendix A 


1. Deriving mean values in the PBL from 'surface data'. 

Prom the earlier version of the GLAS GCM (Tsang and 
Karn, 1973). the 'surface data', u , T and q Q , and inter- 
facial fluxes such as evaporation E, sensible heat flux H 
and friction velocity u # among others are available. With 
Deardorff's parameterization (1972), one can calculate the 
mean values of u^, T m and q m within the PBL by using the 
surface data that are the standard output of a GCM run. 

According to the definition given in chapter 3, 


F w =E/p a (A. 1.1) 

< F.ji=H/(5 > a Cp) (A. 1.2) 

and F=F T +0.6lS m F w (A. 1.3) 

As a first approximation, the S m can be approximated by 

S n =T a ( 1000/P)* 288 (A. 1.4) 

Three stability cases can be determined from the magnitude 
of F, i.e. 

Unstable, if F>0 (A. 1.5) 

Stable, if F<0 (A. 1.6) 

Neutral, if F=0 (A. 1.7) 


If F = 0, the neutral case occurs and the relation be- 
tween the mean atmospheric conditions in the PBL and the 
surface data is given by 
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and 


( V u a>/ u *= 8-4 

(S vnT s va> u ./ I '= 7 -? 
from which, u ffl and S vm can be solved explicitly. 


21 1 
(A. 1.8) 
(A.1 .9) 


If F < 0, the stable case occurs and the equations are 
( V u a )/u * = 8 -4+.6h/L m (A. 1.10) 

(S vm- 2 * * S * va )u * /P= 7.3+. 6h / L m (A. 1.11) 


and 


\= p > 


(A. 1.12) 


from which, u ffl , S ym and can be solved expliciltly. 

If P > 0, the unstable case occur and equations are 

,-.16 


(< V u a )/u * = 8.4(1.-50h/L m )' 

u,/F= 7 . 3( 1 • -5 • 8h/l m )‘ • 47 


and 


vm va 
L m= u * 3s vm /(k « P) 


(A. 1.13) 
(A. 1.14) 
(A.1 .12) 


from which, u m , S Vffl and L m must be solved numerically. 

After having the solution u ffl , S vm and L m , the S m and q m 
can be calculated by the following relations: 

(A. 1.15) 


and 


<V S a>/V( s vm- s va>' P 
(< lm- < la>/ P w = ( S v™-Sv.)/5' 


vm va' 


(A. 1 .16) 


2. Adjustment of net radiation: 


The net radiation from the GCM must be adjusted for the 

GHM to use under the noninteractive mode in order to ac- 

count for the difference in the subgrid parameterization of 

albedo between the value furnished by GCM and the value 

calculated by the GHM formulation. Prom the GCM, one has 
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R *=S *+R t - 4 

n w L gs 


(A. 2.1 ) 


where S y is the absorbed solar radiation in the GCM, R^ is 
the downward long wave radiation and Tg S is the ground sur- 
face temperature provided by the GCM. The downward solar 


radiation S can be calculated from 
w 


V S w7('- a l'> (A. 2. 2) 

♦ 

where a^ is the surface albedo given in the GCM. The net 
radiation in the GEM is given by 




(A. 2.3) 


where T cq is either bare soil surface temperature or canopy 
temperature from the GHM and a^ is the GHM albedo given in 
Section 2.6. 
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